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The maximum efficiency of nano heat engines depends on more than temperature 


Mischa Woods/’^’^ Nelly Ng,^’^ and Stephanie Wehner^’^ 

'University College of London, Department of Physics & Astronomy, London WCIE 6BT, United Kingdom 
^QuTech, Delft University of Technology, Lorentzweg 1, 2611 CJ Delft, Netherlands 
^ Centre for Quantum Technologies, National University of Singapore, 117543 Singapore 

Sadi Carnot’s theorem regarding the maximum efficiency of heat engines is considered to be of fundamental 
importance in thermodynamics. This theorem famously states that the maximum efficiency depends only on the 
temperature of the heat baths used by the engine - but not the specific details on how these baths are actually 
realized. Here, we show that at the nano and quantum scale, this law needs to be revised in the sense that more 
information about the bath other than its temperature is required to decide whether maximum efficiency can 
be achieved. In particular, we derive new fundamental limitations of the efficiency of heat engines that show 
that the Carnot efficiency can only be achieved under special circumstances, and we derive a new maximum 
efficiency for others. This new understanding of thermodynamics has implications for nanoscale engineering 
aiming to construct very small thermal machines. 


Nicolas Leonard Sadi Carnot is often described as the “fa¬ 
ther of thermodynamics”. In his only publication in 1824 [1] 
, Carnot gave the first successful theory of the maximum effi¬ 
ciency of heat engines. It was later used by Rudolf Clausius 
and Lord Kelvin to formalize the second law of thermody¬ 
namics and define the concept of entropy [2, 3] . In 1824 he 
concluded that the maximum efficiency attainable did not de¬ 
pend upon the exact nature of the working fluids [1 ] : 

The motive power of heat is independent of the 
agents employed to realize it; its quantity is fixed 
solely by the temperatures of the bodies between 
which is effected, finally, the transfer of caloric. 

For his “motive power of heat”, we would today say “the effi¬ 
ciency of a reversible heat engine”, and rather than “transfer of 
caloric” we would say “the reversible transfer of heat.” Carnot 
knew intuitively that his engine would have the maximum effi¬ 
ciency, but was unable to state what that efficiency should be. 
Working fluids refers to the substance (normally gas or liquid) 
which is at the hot or cold bath temperatures. 

Carnot also defined a hypothetical heat engine (now known 
as the Carnot engine) which would achieve the maximum ef¬ 
ficiency. Later, this efficiency - now known as the Carnot effi¬ 
ciency (C.E.) - was shown to be 


jyc = 1 - 


/?Hot 
/^Cold ’ 


( 1 ) 


where ficoU, /^Hot are the inverse temperatures of the cold and 
hot baths respectively. 

Unlike the large scale heat engines that inspired thermody¬ 
namics, we are now able to build nanoscale quantum machines 
consisting of a mere handful of particles, prompting many 
efforts to understand quantum thermodynamics (see e.g. [4- 
19]). Such devices are too small to admit statistical methods, 
and many results have shown that the workings of thermody¬ 
namics become more intricate in such regimes [4-7, 15] . 

We show in this report that unlike at the macroscopic 
scale - at which Carnot’s fundamental results hold - there are 
new fundamental limitations to the maximal efficiency at the 
nanoscale. Most significantly, this new efficiency depends on 
the working substance. We find that the C.E. can be achieved. 


but only when the working substance is of a particular form. 
Otherwise, a reduced efficiency is obtained, highlighting the 
significant difference in the performance of heat engines as 
our devices decrease in size. 


Work in the nanoregime 

The basic components of a heat engine (H.E.) are de¬ 
tailed in Eig. 1. The definition of work when dealing with 
nanoscopic quantum systems has seen much attention lately 
[4-9] . Performing work is always understood as chang¬ 
ing the energy of a system, which we call battery. In the 
macroregime, one often pictures raising a weight on a string. 
In the nanoregime, this corresponds to changing the energy of 
a quantum system by lifting it to an excited state (see Eig. 2). 

One aspect of extracting work W is to bring the bat¬ 
tery’s initial state to some final state such that W = 
— tr(p^iTw) > 0. However, a change in energy 
alone, does not yet correspond to performing work. It is im¬ 
plicit in our macroscopic understanding of work that the en¬ 
ergy transfer takes place in an ordered form. When lifting a 
weight, we know its final position and can exploit this precise 
knowledge to transfer all the work onto a third system without 
- in principle - losing any energy in the process. In the quan¬ 
tum regime, such knowledge corresponds to p^ being a pure 
state. When p^ is diagonal in the energy eigenbasis of 
then p^ is an energy eigenstate. We can thus understand work 
as an energy transfer about which we have perfect informa¬ 
tion, while heat, in contrast, is an energy transfer about which 
we hold essentially no information. Clearly, there is also an 
intermediary regime in which we transfer energy, while hav¬ 
ing some - but not perfect - information. 

To illustrate this idea, consider a two-level system battery, 
where we extract work by transiting from an initial energy 
eigenstate \E^){E^\ to another energy eigenstate \E^){E^\, 
where E^ — E^ > 0. Changing the energy, while hav¬ 
ing some amount of information corresponds to changing the 
state of the battery to a mixture p^ = (1 — e)\E^){E^\ + 
£\E^){E^\ for some parameter e G [0,1]. The case of e = 0 
corresponds to doing perfect work. The smaller e is, the closer 
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FIG. 1. A heat engine extracts work from the temperature differ¬ 
ence of a hot bath (red) at temperature Thoi and cold bath (blue) at 
temperature Tcoid- The term bath indicates that the initial states of 
these systems are thermal, rnot = exp(—/3Hot(THot)/^Hot and rcoid = 
exp(—/lcoidTfcoid)/^Coid with inverse temperatures /(not = l/Tkot 
and /3coid = l/Tcoid, and partition functions Zuot and Zcoid respec¬ 
tively. The machine itself corresponds to a quantum system M with 
Hamiltonian Hm, starting in an arbitrary state The battery in¬ 
dicates the system in which we will store work and is illustrated in 
detail (Figure 2). Let be the starting state of the work system, 
and Hv/ its Hamiltonian. Operating the heat engine for one cycle for 
some particular time t, corresponds to applying a unitary transform 
U(t) to the both baths, the actual machine, and also the work sys¬ 
tem. In order to account for all energy transfers we will demand that 
[(7(f), H] = 0, where H = Hu + TTcoid + fTnot + ^w. That is, 
U (f) conserves total energy. We will furthermore demand that op¬ 
erating the heat engine does not affect the actual machine, i.e., the 
final state of the machine pjj = p^, but the choice of p^ may af¬ 
fect the efficiencies we can obtain. In the language of [5] , we thus 
treat the machine as a catalyst. The efficiency rj = lLext/A77 of 
the machine measures the amount of work extracted H4xt in relation 
to the energy invested as measured by the energy change of the hot 
bath Ai7. Given two temperatures Tnot and Tcoid, finding the opti¬ 
mal heat engine for such temperatures would correspond to deciding 
on the best machine (p^, Hu), evolution U (f), and bath structures 
Huot and Hcou- To understand the role of the structure of the cold 
bath, we here fix T/coid while allowing an otherwise arbitrary choice 
of Hnot- 


we are to the situation of perfect work. One can characterize 
this intermediary regime by the von Neumann entropy S (p^)- 
For perfect work, S{p\i) = 0, while for heat transfer, under 
a fixed average energy, the two-level battery becomes ther¬ 
mal, since the thermal state maximizes entropy for a fixed 
energy [20] . Intuitively, we may also think about e as the 
failure probability of extracting work. When e > 0, what is 
relevant is not e as an absolute, but relative to the energy Wgxt 
that is extracted. We are thus interested in AS'/Wext where 
AS” = S'(pw) ~ '^(Pw) *^he change in entropy of the battery. 
For our investigations the limit AS'/Wext —f 0 will be of par¬ 
ticular interest, and corresponds to performing near perfect 
work. Our analysis applies to arbitrarily small heat engines, 



FIG. 2. A battery is a work-storage component of a heat engine. 
In the nanoregime, a minimal way of modeling the battery is as a 
two-level system [5] . Performing work corresponds to “lifting” the 
state from ground state to the excited state, where the energy gap is 
fine-tuned to the amount of work Wext to be done. While an arbi¬ 
trary energy spacing is difficult to realize in a two-level system, it 
can be done by picking two levels with the desired spacing from a 
quasi-continuum battery; this battery comprises of a large but finite 
number of discrete levels which form a quasi-continuum. Such a bat¬ 
tery closely resembles the classical notion of a “weight attached to 
a string” as considered in [21] . The battery can be charged from a 
particular state (e.g. the ground state) to any of the higher levels. 


even if this machine was run for only one cycle. We empha¬ 
size that this is not a restriction of the analysis, but rather a 
strong and appealing feature because it is indeed the relevant 
case when we consider few qubit devices, and a small number 
of experimental trials. 


No perfect work 

Before establishing our main result, we first show that in 
the nanoscopic regime, no heat engine can perform perfect 
work (e = 0). That is, the efficiency of any such heat engine 
is zero. More formally, it means that there exists no global 
energy preserving unitary (see Fig. 1) for which W/xt > 0. 

Efficiency 

Clearly, however, heat engines can be built, prompting the 
question how this might be possible. We show that for any 
e > 0, there exists a process such that Wext > 0. There¬ 
fore, a heat engine is possible if we ask only for near perfect 
work. Interestingly, even in the macroscopic regime, we can 
envision a heat engine that only extracts work with probability 
1 — e, but over many cycles of the engine we do not notice this 
feature when looking at the average work gained in each run. 

To study the efficiency in the nanoscale regime, we make 
crucial use of the second laws of quantum thermodynamics [5] 
. It is apparent from these laws that we might only discover 
further limitations to the efficiency than we see at the macro¬ 
scopic scale. Indeed they do arise, as we find that the effi¬ 
ciency no longer depends on just the temperatures of the heat 
baths. Instead, the explicit structure of the cold bath Hamil¬ 
tonian F/coid becomes important (a similar argument can be 
made for the hot bath). Consider a cold bath comprised of n 
two-level systems each with its own energy gap, where n can 
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be arbitrarily large, but finite. Let us denote the spectral gap 
of the cold bath—the energy gap between its ground state and 
first excited state—by Smin- We can then define the quantity 

.2'> 

—j— g /3coldt?niin ' 

and study the efficiency in the quasi-static limit. This means 
that the final state of the cold bath is thermal, and its final tem¬ 
perature Tf is higher than Tgoid by only a positive infinitesimal 
amount. 

Whenever < 1, we show that the maximum and attain¬ 
able efficiency is indeed the familiar C.E., which can be ex¬ 
pressed as 


V = 



/?Hot \ 

/^Coid — Pmt J 


-1 


( 3 ) 


However, when H > 1, we find a new nanoscale limitation. In 
this situation, the efficiency is only 


?7 = 



PCold — PHot / 


-1 


( 4 ) 


for a quasi-static heat engine. One might hope to obtain a 
higher efficiency compared to (4) by going away from the 
quasi-static setting, however we also show that such an effi¬ 
ciency is always strictly less than the C.E. 

The restriction of near perfect work per cycle can now be 
further justified by examining how well the heat engine per¬ 
forms when the machine runs over many cycles: we find that 
if H < 1, the heat engine can be run quasi-statically with an 
efficiency arbitrarily close to the C.E. while extracting any fi¬ 
nite amount of work with an arbitrarily small entropy increase 
in the battery. 


Comparison with standard entropy results 


Eor any system in thermal contact with a bath at tempera¬ 
ture T, consider the Helmholtz free energy F{p) — Tr{Hp) — 
kTS{p), where S{p) = —tr(plnp) is the von Neumann en¬ 
tropy of p. In the macroregime, the usual second law states 
that the Helmholtz free energy never increases. 
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FIG. 3. For fixed temperatures Tcom, Tuot, the efficiency of a 
nanoscale heat engine depends on the structure of the cold hath. 

At the nano/quantum scale, Carnot’s statement about the universality 
of heat engines does not hold. We find that the maximum efficiency 
of a heat engine, does depend on the “working fluid”. In (a) the en¬ 
ergy gaps are small enough to allow the heat engine to achieve C.E., 
i.e., H < 1. In (b) the efficiency of the heat engine is reduced below 
the C.E. because the energy gap of the qubits are above the critical 
value H > 1. 


catalysis [5, 22], however, this corresponds to a degradation 
of the machine in each cycle. 

It is illustrative to analyze our problem when we apply just 
the standard second law in (5) to derive bounds on the effi¬ 
ciency, which is indeed a matter of textbook thermodynam¬ 
ics [23] . However, we here apply the law precisely to the heat 
engine model as given in Eig. 1, in which all energy flows 
are accounted for and (near) perfect work is performed. One 
might wonder whether the limitations we observe are just due 
to an inaccurate model or our demand for near perfect instead 
of average work, and might thus also arise in the macroregime. 
That is, are these newfound limitations really a consequence 
of the need to obey a wider family of second laws, or would 
the standard free energy predict the same things when the en¬ 
ergy is quantized? 


F{po) > F{p,) , (5) 

when the system goes from a state po to a state pi. This, 
however, is but one of many conditions necessary for a state 
transformation [5] . The limitations we observe are a conse¬ 
quence of the fact that in the nanoregime, possible transitions 
are governed by a family of generalized second laws. The fact 
that more laws appear in this regime can intuitively be under¬ 
stood as being analogous to the fact that when performing a 
probabilistic experiment only a handful of times, not just the 
average, but other moments of a distribution become relevant. 
Indeed, all second laws converge to the standard second law 
in the limit of infinitely many particles [5], illustrating why 
we are traditionally accustomed to only this second law. The 
standard second law also emerges in some regimes of inexact 


We show independently of whether we consider perfect 
(e = 0) or near perfect (AS'/Wext —> 0) work - that according 
to the standard free energy in (5), the maximum achievable ef¬ 
ficiency is the C.E. Eurthermore, we recover Carnot’s famous 
statement that the C.E. can be achieved for any cold bath (i.e. 
for a cold bath with any finite dimensional pure point spec¬ 
trum). We also see that C.E. is only achieved for quasi-static 
H.E.s. We prove this without invoking any additional assump¬ 
tions than those laid out here, such as reversibility or that the 
system is in thermodynamic equilibrium at all times. There¬ 
fore, with our setup we recover exactly what Carnot predicted, 
namely that the maximum efficiency of the H.E. is indepen¬ 
dent of the working substance. This rules out that our inabil¬ 
ity to achieve what Carnot predicted according to the macro¬ 
scopic laws of thermodynamics is not only a consequence of 
an overly stringent model. 
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Extensions to the setup 

One could ask whether at the nano regime, if a less strin¬ 
gent model would allow one to recover Carnot’s predictions. 
Specifically, what if we consider any final state of the bat¬ 
tery which is e away in trace distance from the desired 
final battery state In this case, we show that as 

long as one still considers the extraction of near perfect work 
AS'/Wext —> 0, our findings remain unchanged: when O > 1, 
C.E. cannot be obtained. 

In a similar vein, one could imagine that the final compo¬ 
nents of the heat engine become correlated between them¬ 
selves, and that this would allow one to always achieve the 
C.E.. According to macroscopic laws of thermodynamics, 
correlations between the final components always inhibit one 
from achieving the C.E.. We show that at the nanoscale such 
correlations can also be ruled out as a means to achieve the 
C.E. when 17 > 1. 

These results thus show the inevitability that the maximum 
efficiency of a nanoscale heat engine depends on more infor¬ 
mation about the thermal baths rather than just the tempera¬ 
ture. 
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Conclusion 

Our work establishes a fundamental result for the opera¬ 
tion of nanoscale heat engines. We find all cold baths can 
be used in heat engines. However, for all temperatures Tcoid 
and Thoi of the cold and hot baths, there exists an energy gap 
£'min(/3coidi i^Hot) of the two-level systems forming the cold 
bath above which the optimal efficiency is reduced below the 
C.E.. Viewed from another direction, for a fixed energy gap 
-E'niin(/3coid)/^Hot), whether the C.E. can be achieved depends 
on the relation between Tkot and Tcoid as illustrated in Eig. 4. 
Loosely speaking, the C.E. can be achieved whenever the two 
temperatures are unequal but not too far apart. One might 
wonder why this restriction has not been observed before in 
the classical scenario. There, the energy spectrum is continu¬ 
ous or forms a quasi-continuum, and hence we can only access 
the C.E. regime. 

Our result is a consequence of the fact that the second law 
takes on a more complicated form in the nanoregime. Next 
to the standard second law, many other laws become rele¬ 
vant and lead to additional restrictions. Erom a statistical 
perspective, small numbers require more refined descriptions 
than provided by averages, and as a result thermodynamics 
becomes more complicated when considering systems com¬ 
prised of few particles. Similar effects can also be observed in 
information theory, where averaged quantities as given by the 
Shannon entropy need to be supplemented with refined quan¬ 
tities when we consider finitely many channel uses. 

In the macroscopic regime, for completeness, we ruled out 
the possibility that the observed limitations on efficiency is 
a consequence of our demand for near perfect work, or the 
fact that we are using systems with discrete (sufficiently large 
spaced) spectra. This verification was achieved by showing 
that the C.E. can indeed always be attained (regardless to the 
size of an energy gap if present) when extracting near per- 


FIG. 4. Comparison of the nano/quantum-scale efficiency versus 
macroscopic efficiency (C.E.), in the quasi-static regime. Top: 

Efficiency vs. the energy gap E^m of ffcoid- According to Eqs. (3), 
(4), for any /^Hot < Pcou (Tcoid < Tkot) one can achieve C.E. when 
the structure JTcoid of the cold bath has sufficiently small energy gap. 
If the energy gap is too large, we find a reduced efficiency. This has 
been plotted for Tnot = 15 and Tcoid = 10. Middle: Efficiency vs. 
Tcoid. For every Hcou, there exists a temperature regime (Tcow vs. 
Thoi) such that C.E. cannot be achieved. This happens as Tcoid gets 
further from the temperature of the hot bath Thm = 20, where we 
have used E,„in = 15. Bottom: Efficiency vs. Tnot. Similarly, we see 
this feature of not being able to achieve C.E., as the temperature of 
the hot bath increases relative to Tcoid = 5, where again Emin ~ 15. 


feet work, when we are in such large systems that only the 
standard second law is relevant. One might wonder whether 
heat engines that do not operate quasi-statically, or employ¬ 
ing quantum coherences would allow us to achieve the C.E. 
independent of the structure of the cold bath. As we show in 
the Supplementary Material, both do not help. 

There are several works [15, 21, 24-26] that have analyzed 
the efficiencies of heat engines and obtained C.E. as the max¬ 
imal efficiency. Common to all these approaches is that they 
consider an average notion of work, without directly account¬ 
ing for a contribution from disordered energy (heat). Instead, 
one keeps the entropy of the battery low [21] , or bound the 
higher moments of the energy distribution [25] . These only 
limit contributions from heat, but do not fully prevent them. 
Our notion of (near) perfect work now makes this aspect of 
macroscopic work explicit in the nanoregime. Needless to 
say, imperfect work with some contribution of heat can also 
be useful. Yet, it does not quite constitute work if we cannot 
explicitly single out a contribution from heat. One could con¬ 
struct a machine which extracts some amount of energy, with 
some non-negligible amount of information. We can prove in 
this case that Carnot’s efficiency can even be exceeded [27] . 
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This should not come as a surprise, because we are no longer 
asking for work - energy transfer about which we have (near) 
perfect information. 

Our work raises many open questions. We see that the 
quasi-static efficiency has a discontinuous derivative with re¬ 
spect to Thoi, Tcoid or E'niin at O = 1, as illustrated in Fig. 4, 
which is often associated with a phase transition. It is unclear 
whether this phenomenon can also be understood as a phase 
transition - absent at the macroscopic scale - and whether there 
is an abrupt change in the nature of the machine when crossing 
the 0 = 1 boundary. 

It would furthermore be satisfying to derive the explicit 
form of a hot bath, and machine attaining Carnot - or new 
Carnot - efficiency. One might wonder whether a non-trivial 
machine is needed at all in this case. To illustrate 

the dependence on the bath, it was sufficient to consider a bath 
comprising solely of qubits. The tools proposed in the Supple¬ 


mentary Material can also be used to study other forms of bath 
structures, yet it is a non-trivial question to derive efficiencies 
for such cold baths. 

Most interestingly, there is the extremely challenging ques¬ 
tion of deriving a statement that is analogous to the C.E., but 
which makes explicit the trade-off between information and 
energy for all possible starting situations. In a heat engine, we 
obtain energy from two thermal baths about which we have 
minimal information. It is clear that the C.E. is thus a spe¬ 
cial case of a more general statement in which we start with 
two systems of a certain energy about which we may have 
some information, and we want to extract work by combining 
them. Indeed, the form that such a general statement should 
take is by itself a beautiful conceptual challenge, since what 
we understand as efficiency may not only be a matter of work 
obtained vs. energy wasted. Instead, we may want to take a 
loss of information about the initial states into account when 
formulating such a fully general efficiency. 
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Methods 


Setup of heat engine The workings of a heat engine have been 
described in Fig. 1, which we expand in mathematical detail 
here. Consider the initial global system 

PColdHotMW = Pcold ® Phoi Pm ® Pw- (6) 

The hot bath Hamiltonian Hhoi can be chosen arbitrarily, as 
long as p^Q( is the corresponding thermal state of temperature 
Thoi- Similarly, the machine and its Hamiltonian (p^,FfM) 
can be chosen arbitrarily. Given any cold bath (Pco^, fTcoid) 
such that Pcoi(j is a thermal state at temperature Tcoid < Tnot, 
one can extract work and store it in system W. This process 
then corresponds to inducing the transition 

PColdHotMW ^ /^ColdHotMW J (^) 

where trHot(PcoidHotMw) = PLidw ® Pm- We have p^ = pj^, 
that is, the machine is not degraded in the process. This also 
means that we preserve the tensor product structure between 
ColdW and M: if the machine is initially correlated with some 
other system, we do not want to destroy such correlations 
since this would also mean that we degrade the machine. 

To quantify the amount of extractable work, we apply the 
generalized second laws derived in [5] . The initial cold bath 
Pcoid thermal, and therefore diagonal in the energy eigen- 
basis, while the initial battery state p^ is also a pure energy 
eigenstate (see Fig. 2). Since here energy conserving U{t) 
never increase coherences between energy eigenstates [5] , 
we can therefore conclude that p^oidw diagonal in the 

energy eigenbasis. We can thus invoke the necessary and 
sufficient conditions for a transformation to be possible [5] . 
Specifically, pji^y (g) p^ ^ p^„yw iff Va > 0, 

(PCold ® PW J '^ColdW ) — (PColdW J '^ColdW ) > (^) 

where is the thermal state of the joint system (cold bath 

and battery) at temperature T^ot- The generalized free energy 
Fa is defined as 


Pa{p,T) 


1 

PhoI 


[D„(p||r) - InZHot] , 


(9) 


where Z?Q,(p||r) are known as a-Renyi divergences. For states 
p, T which are diagonal in the same eigenbasis, the Renyi di¬ 
vergences can be simplified to 

DaipWr) = -^—r ln^p“gi-“, (10) 

a — 1 


where pi, qi are the eigenvalues of p and r respectively. The 
case a = 1 is defined by continuity in a. Taking the limit 
a —> 1 for (9), one recovers the Helmholtz free energy, 
F(p) = (fT)p — /^Hol‘^(p)- TFsing the second laws [5] is 
a powerful tool, since when searching for the optimum effi¬ 
ciency, we do not have to optimize explicitly over the possible 
machines (pM, Hm), the form of the hot bath Haot, or the en¬ 


ergy conserving unitary U{t). Whenever (8) is satisfied, then 
we are guaranteed a suitable choice exists and hence we can 
focus solely on the possible final states Pcoidw- 

Since we know that Pcoidw ts a diagonal in the energy 
eigenbasis, the correlations between cold bath and battery can 
only be classical (w.r.t. energy eigenbasis). However, even 
such correlations cannot improve the efficiency; we show in 
the Supplementary Material that we may take the output state 
to have the form Pcoidw = Pcoid ® Pw order to achieve 
the maximum efficiency. According to Fig. 2, consider p^ = 

I-^wXT'wI Pw = (1 ~ ‘S)|T'w)(-®wl "f ^iT'wXf^wl- From 
the second laws (8), we see that the maximum amount of 
extractable work is given by the largest value of IFext = 
such that the state transition Pcoid®Pw Pcoid®Pw 
is possible. The form of IFext (derived in the Supplementary 
Material) is 

kFext = inf Wa, 
a>0 

Wa = ^^ -^[ln(A-£“)-Q!ln(l-e)], 

PHot(Q! — 1) 

A = _ 

E /fy 1 —a ’ 

ip [ 

where pi = ^ m ^ ^ are probabilities of the 

^ ^/3cold ^ ^/3 ho. ^ 

thermal state of the cold bath at temperatures Tcoid,THot re¬ 
spectively, and p' are the probability amplitudes of state Pcoid 
when written in the energy eigenbasis of HcoM- The quantity 
Wext is dependent on the initial and final cold bath Pcojd, Pcoid’ 
the hot bath temperature Thoi, and the allowed failure proba¬ 
bility e. The main difficulty of evaluating Wgxt comes from 
the infimum over a, - indeed we know examples in which it 
can be obtained at any a > 0 depending on Paot, ,5coid and 
other parameters. 


( 11 ) 

( 12 ) 

(13) 


The efficiency p, however, is not determined the maximum 
extractable work, but rather by a tradeoff between IText and 
the energy drawn from the hot bath. More precisely. 


p := 


M4, 


AHot’ 


(14) 


where 


AHot := tr(fTHotpSot) “ ti'(-^Hotpiot) (15) 


is the mean energy drawn from the hot bath. Since 
-ffcoidHotMW = ^Coid + -ffHot + iTM + -ffw is void of interaction 
terms, and since total energy is preserved, we can also write 
the change of energy in the hot bath, in terms of the energy 
change in the remaining systems. That is. 


AHot = ACold -f AfF. (16) 


where 


ACold := tr [iTcoidPcoid] “ tr [-^CoidPcoid] > (17) 
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and 

AW" := ti-(i7wPw) - tr(i7w/9^). (18) 

is the change in average energy of the cold bath and battery. 
We thus see that the efficiency can be described wholely in 
terms of the battery and the cold bath. 

Macroscopic second law We first analyze the efficiency in the 
macroscopic regime, where only the usual free energy (a = 1) 
dictates if a certain state transition is possible. The main ques¬ 
tion is then: given an initial cold bath Hamiltonian 7/coid, what 
is the maximum attainable efficiency considering all possible 
final states Pcoid^ cases of perfect and near perfect 

work, we find that the efficiency is only maximized whenever 
Pcoid 1® (^) ^ thermal state, and (B) has a temperature /?/ arbi¬ 
trarily close to /3coid- We refer to this situation as a quasi-static 
heat engine. Moreover, we find that the maximum is the C.E. 
and that it can be achieved by any given iTcoid- These results 
rigorously prove Carnot’s findings when only the usual free 
energy is relevant. 

Nanoscale second laws Here, when considering perfect work, 
we are immediately faced with an obstacle: the constraint at 
a = 0 implies that Wext > 0 is not possible, whenever 
is of full rank. This is due to the discontinuity of Dq in the 
state probability amplitudes, and is similar to effects observed 
in information theory in lossy vs. lossless compression: no 
compression is possible if no error however small is allowed. 
However, when considering near perfect work, the Dq con¬ 
straint is satisfied automatically. We thus consider the limit 
AS/kkext ^ 0. 

The results for the macroscopic second law implies an up¬ 
per bound on both the maximum extractable work and ef¬ 
ficiency for nanoscopic second laws, since the constraint of 
generalized free energy at a = 1 is simply one of the many 
constraints described by all a > 0. It thus follows from the 
macroscopic second law results, that if we can achieve the 
C.E., we can only do so when both (A) and (B) are satisfied. 
Consequently, we analyze the quasi-static regime. Eurther- 
more, we specialize to the case where the cold bath consists 
of multiple identical two-level systems, each of which are de¬ 
scribed by a Hamiltonian with energy gap E. 

Eirstly, we identify characteristics that £ should have, such 
that near perfect work is extracted in the limit Pf —t /3coid (i-e. 
when (A) and (B) are satisfied). We then show two technical 
results: 

1. The choice of e (as a function of Pf) simplifies the min¬ 
imization problem in (11), by reducing the range the 
variable a appearing in the optimization of Wext- Un¬ 
der the consideration of near perfect work, e can be cho¬ 
sen such that the optimization of a is over a > k for 
some K G (0,1], instead of a > 0. The larger k is for a 
chosen e, the slower AS/Wgxt converges to zero. 

2. We analyze the following cases separately: 

• Eor H < 1, e can always be chosen such that the 
infimum in (11) is obtained in the limit a 1. 
Evaluating the efficiency in the limit a —> 1 cor¬ 
responds to the C.E.. 


• Eor fl > 1, we show that for the best choice 
of s, the infimum in (11) for W"ext is obtained at 
a —7^ oo. Eurthermore, fl > 1 means that up to 
leading order terms, Wi > Woo for Wa defined 
in (12). But we know that the quantity Wi gives 
us C.E.. Therefore, the efficiency is strictly less 
than Carnot. 
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In this Supplementary Material, we detail our hndings. Sections A-C are aimed at giving the reader an overview of the 
important concepts regarding heat engines, and to introduce the quantities of interest. Firstly, in Section A we describe the setup 
of our heat engine, the systems involved, and how work is extracted and stored. By using this general setup, we then proceed in 
Section B to introduce conditions for thermodynamical state transitions in a cycle of a heat engine. In Section C, we introduce 
the formal dehnition of efficiency, and specify how can this quantity be maximized over a set of free parameters (involving the 
bath Hamiltonian structure). 

After providing these guidelines, we start in Section D to apply the macroscopic law of thermodynamics. We have performed 
the analysis with the generalization of allowing for an arbitrarily small probability of failure. The results in this section might 
be familiar and known to the reader, however from a technical perspective, their establishment is helpful for proving our main 
results (in Section E) about nanoscale systems. In Section E, we apply the recently discovered generalizations of the second law 
for small quantum systems. The results in Section D and Section E are summarized at the beginning of each section, for the 
reader to have a concise overview of the distinction between thermodynamics of macroscopic and nanoscopic systems. Finally, 
in Section F, we show that even when considering a more general setup, these results obtained in Section E remain unchanged. 
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A. The general setting for a heat engine 


A heat engine is a procedure for extracting work from a temperature difference. It comprises of four basic elements: two 
thermal baths at distinct temperatures Thoi and Tcoid respectively, a machine, and a battery. The machine interacts with these 
baths in such a way that utilizes the temperature difference between the two baths to perform work extraction. The extracted 
work can then be transferred and stored in the battery, while the machine returns to its original state. 

In this section, we describe a fully general setup, where all involved systems and changes in energy are accounted for explicitly. 
Let us begin with the total Hamiltonian 


Ht — Tfcold + Hy{oi + Hyi + Tfwj 


(Al) 


where the indices Hot, Cold, M, W represent a hot thermal bath (Hot), a cold thermal bath (Cold), a machine (M), and a battery 
(W) respectively. Let us also consider an initial state 

PColdHotMW = '^Cold ® '^Hot ® Pm® P%- (A2) 


The state (''"coid) '^^e initial thermal state at temperature Thoi (Tcoid), corresponding to the hot (cold) bath Hamil¬ 

tonians HRotjHcoid- More generally, given any Hamiltonian H and temperature T, the thermal state is defined as t = 
/ksT^ For notational convenience, we shall often use inverse temperatures, defined as /3h ■= l/fcsTHot and 
Pc '■= l/fcsTcoid where Ub is the Boltzmann constant. Without loss of generality we set Tcoid < Thoi- The initial machine 
(p^, Hm) can be chosen arbitrarily, as long as its final state is preserved (and therefore the machine acts like a catalyst). 

We adopt a resource theory approach, which allows all energy-preserving unitaries U{t) on the global system, i.e. such 
unitaries obey \U{t), ILcoidHotMw] = 0. If ^Hot) and (p^, Hm) can be arbitrarily chosen, these correspond to the set of 
catalytic thermal operations [."i] one can perform on the joint state ColdW. This implies that the cold bath is used as a resource 
state. By catalytic thermal operations that act on the cold bath, using the hot bath as a thermal reservoir, and the machine as a 
catalyst, one can possibly extract work and store it in the battery. 

The aim is to achieve a final reduced state PcoidHotMW’ 


PColdMW — THot(PcoldHotMw) — PCoid ® Pm® Pw^ 


(A3) 


where pj^ = p^, i.e. the machine is preserved, and Pcoid? Pw '^^e final states of the cold bath and battery. In Section F, we 
will consider the case in which there are correlations between the final state of the cold bath, hot bath, battery and or machine. 
We will find that the correlations do not change our results. For any bipartite state pab, we use the notation of reduced states 
PA := trB(pAB), Pb := trA(pBA)- 

Finally, we describe the battery such that the state transformation from PcoidHotMW PcoidHotMW stores work in the battery. 
This is done as follows: consider the battery which has a Hamiltonian (written in its diagonal form) 




(A4) 


For some parameter e G [0,1), we consider the initial and final states of the battery to be 

p°^=\E,){E,\^ (A5) 

Pw = (1 ~ £)\Ek){Ek\w + £\Ej){Ej\^ff (A6) 

respectively. The parameter Wext is defined as the energy difference 

VFext := E^ - Ej. (A7) 

where we define E'^ > E^ such that ILext > 0. In the case where Wext is a value such that the transition Pc^idw Pcoidw 
is possible via catalytic thermal operations, it corresponds to extracting work. We refer to the parameter e as the probability of 
failure of work extraction. Note that e in Eq. (A6) is also the trace distance 

^(PW) = ^l|p-cr||i (A8) 

between \Ej){Ej\^ and \Ek){Ek\^f^■ In Section F, we will generalize this definition to include all final states of the battery p^. 
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FIG. 5. The setting of a working heat engine. 


which are a trace distance e from the ideal final battery state \Ek){Ek\^. We show that our findings regarding the achievability 
of C.E. remains unchanged. 

Throughout our analysis, we deal with two distinct scenarios of work extraction as defined below. 

Definition 1. (Perfect work) An amount of work extracted Wext is referred to as perfect work when e = 0. 

The next definition of work involves a condition regarding the von Neumann entropy of the final battery state. Let AS' be the 
von Neumann entropy of the final battery state. When the initial state is pure, we have 

AS := -tr(pwlnpw)- (A9) 

When the final battery state is given by Eq. (A6), its probability distribution has its support on a two-dimensional subspace of 
the battery system, this definition also coincides with the binary entropy of £, 

h 2 (£) = —eln£ — (1 — e) ln(l — e) = AS. (AlO) 


Definition 2. (Near perfect work) An amount of work extracted Wext A referred to as near perfect work when 


1) 0 < e < I, for some fixed I < 1 and 


2 ) 0 < 


AS 

W^x, 


< p for any p > 0, i.e. 


AS 


is arbitrarily small. 


In the main text, we have provided a detailed discussion regarding the physical meaning of perfect work and near perfect 
work, and the necessity for considering these quantities. As we will see later in the proof to Lemma 5, 1) and 2) in Def. 2 are 
both satisfied if and only if 


, AS 

lim —— = 0. 
£->0+ Wext 


(All) 


Since the initial state PcoidHotMW A diagonal in the energy eigenbasis, and since catalytic thermal operations do not create 
coherences between energy eigenstates, therefore PcoWMW has to be diagonal in the energy eigenbasis. Furthermore, (as already 
stated above) in Section F, we extend the setup to include correlation in the final state between the battery, cold bath and machine 
and more general final battery states. 

Note that in our model we allow the battery to have arbitrarily many (but finite) eigenvalues. One can compare this to the two- 
dimensional battery used in [5], referred to as the wit. Having a minimal dimension, the wit is a conceptually very useful tool 
to visualize work extraction. However, it has the disadvantage that the energy spacing, i.e. the amount of work to be extracted, 
has to be known a priori to the work being extracted in order to tune the energy gap of the wit. The more general battery, which 
we describe in Eq. (A4), requires a higher system dimension, but has the advantage that it can form a quasi-continuum and thus 
effectively any amount of work (i.e. any Wext > 0) can be stored in it without prior knowledge of the work extraction process. 
We will see that our results are independent of nw > 2. 

To summarize, so far we have made the following minimal assumptions: 

(A.l) Product state: There are no initial nor final correlations between the cold bath, machine and battery. Initial correlations 
we assume do not exist, since each of the initial systems are brought independently into the process. This is an advantage 
of our setup, since if one assumed initial coherence, one would then have to use unknown resources to generate them in 
the first place. We later also show that correlations between the final cold bath and battery do not provide improvements 
in maximum extractable work or efficiency. 
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(A.2) Perfect cyclicity: The machine undergoes a cyclic process, i.e. Pm = Pm- 

(A.3) Isolated quantum system: The heat engine as a whole, is isolated from and does not interact with the world. This assump¬ 
tion ensures that all possible resources in a work extraction process has been accounted for. 

(A.4) Finite dimension: The Hilbert space associated with PcoidHotMW finite dimensional but can be arbitrarily large. Moreover, 
the Hamiltonians f^^coid, ^Hot, Hm and iTw all have bounded pure point spectra, meaning that these Hamiltonians have 
eigenvalues which are bounded. 


After defining the set of allowed operations, and describing the desired state transformation process, one can then ask: what 
conditions should be fulfilled such that there exists a hot bath HHot), and a machine (pm) Hm) such that Pc^yw —>■ Pcoidw 
is possible? Throughout this document we use to denote a state transition via catalytic thermal operations. 

In Section D, by assuming the macroscopic law of thermodynamics governs the heat engine, we derive the efficiency of a 
heat engine, and verify the long known Carnot efficiency as the optimal efficiency. We do this for both cases where £ = 0 and 
when £ is arbitrarily small. In Section E, we analyze the same problem under recently derived second laws, which hold for small 
quantum systems. We show that these new second laws lead to fundamental differences to the efficiency of a heat engine. 

Throughout our analysis, a particular notion that describes thermodynamical transitions will be important towards achieving 
maximum efficiency. We therefore define this technical term, which will be used throughout the manuscript. 


Definition 3. (Quasi-static) A heat engine is quasi-static if the final state of the cold bath is a thermal state and its inverse 
temperature (if only differs infinitesimally from the initial cold bath temperature, i.e. Pf = (3c — g, where 0 < p <?; 1 . 

Since throughout this analysis we frequently deal with arbitrarily small paramaters £, g, we also introduce beforehand the 
notation of order function 0 (a;), o{x), which denotes the growth of a function. 


Definition 4. (Big Q, small o notation [28]) Consider two real-valued functions P{x),Q{x). We say that 
1. P(x) = Q(Q(x)) in the limit x ^ a iff there exists ci, C 2 > 0 and (5 > 0 such that for all lx — a| < 5, c\ < 


2. P(x) = o{Q{x)) in the limit x —> a iff there exists C 3 > 0 such that lim 


P(x) 


Q(x) 


= C3. 


EM 

Qix) 


< C2. 


Remark 1. In Def.4, if the limit of x is unspecified, by default we take a = 0. In [28], these order terms were only defined for 
X ^ 00. However, choosing a general limit x —>■ a can be done by simply defining the variable x' = l/(x — a), and x —> a“'" is 
the same as taking x' —>■ 00. 

We also list a few properties of these functions here for x —0, which will help us throughout the proof: 

a) For any c ^ 0, 0(c • P(x)) = 0(P(x)). 

b) For any functions Pi(x) and P 2 (x), 0(Pi(x)) -h 0(P2(x)) = 0 (max{|Pi(x)|, |P 2 (x)|}). 

c) For any functions Pi(x) and P 2 (x), 0(Pi(x)) ■ 0(P2(^)) = 0{Piix)P2(x)). 

d) For any functions Pi(x) and P 2 {x), 0{Pi{x))/0{P2(x)) = 0{Pi(x)/P 2 {x)). 


Definition 3 has two direct implications for a quasi-static heat engine: 

(i) The temperature of the final state of the cold bath Tf, only increases w.r.t. its initial temperature by an infinitesimal 
amount, i.e. Tf = Tcoid + ^coid 9 + ©(p^)- 

(ii) The amount of work extracted is infinitesimal: as we shall see later, the extractable perfect and near perfect work IFext > 0 

(see Defs. (1), (2)) is of order 0{g). This follows from using Eq. (D5) for the case where ^ thermal state with 

inverse temperature /3/ = /3c — g, and calculating the Taylor expansion of Wext about g = 0. 


B. The conditions for thermodynamical state transitions 


In this section, we briefly state the laws which govern the transitions from initial, PcddHoiMW final, PcoidHotMW states for one 
cycle of our heat engine. By applying these laws, the amount of extractable work W^xt can be quantified and expressed as a 
function of the cold bath. 
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1. Second law for macroscopic systems 


The cold bath, machine and battery form a closed but not isolated thermodynamic system. This means only heat exchange 
(and not mass exchange) occurs between these systems and the hot bath. Therefore, a transition from PcoWMW PcoidMW "'ill 
be possible if and only if the Helmholtz free energy, F does not increase 

^(PColdMw) ^ -^(PColdMw)) (1^1) 

where 

F{p)-.= {H),-^S{p), (B2) 

and S{p) := —tr(plnp) and {H)p := ti{Hp) being the entropy and the mean energy of state p respectively. Throughout the 
manuscript, whenever the state is a thermal state at temperature /3, we use the shorthand notation (fTcoid)/? and S{0). 

The Helmholtz free energy bears a close relation to the relative entropy, 

Zl(p||cr) = tr(plnp) — tr(plncr). (B3) 

Whenever p and a are diagonal in the same basis, the relative entropy can be written as 

Dip\\cr) =^Piln—, (B4) 

where Pi,qi are the eigenvalues of p and a respectively. Now, for any Hamiltonian H, consider rp = jZp, which is the 
thermal state at some inverse temperature /3, with partition function Zp = tr[e“^^], and denote its eigenvalues as qi. Then for 
any diagonal state p with eigenvalues pi, and denoting {Ei}i as the eigenvalues of H, 

D{p\\Tp) = ^piln— = -S{p) F^Pr{l3Ei FlnZp) = l3E{p) FlnZp. (B5) 

Qi 


This implies that 

Pi.P) = ^[D{.p\\tp) - \nZp]. (B6) 

In Section D we will solve Eq. (Bl) in order to evaluate the maximum efficiency. 


2. Second laws for nanoscopic systems 


In the microscopic quantum regime, where only a few quantum particles are involved, it has been shown that macroscopic 
thermodynamics is not a complete description of thermodynamical transitions. More precisely, not only the Helmholtz free 
energy, but a whole other family of generalized free energies have to decrease during a state transition [5]. This places further 
constraints on whether a particular transition is allowed. In particular, these laws also give necessary and sufficient conditions, 
when a system with initial state Pcoidw transformed to final state p^oidw (both diagonal in the energy eigenbasis), with the 

help of any catalyst/machine which is returned to its initial state after the process. 

We can apply these second laws to our scenario by associating the catalyst with and considering the state transition 
Pw ® '’’Cold Pw ® Pcoid described in Section A. Note that the initial state is block-diagonal in the energy 

eigenbasis (for the battery by our choice, and for the cold bath because it is a thermal state). By catalytic thermal operations, the 
final state is also block-diagonal in the energy eigenbasis. Furthermore, according to the second laws in [5], the transition from 
Pw ® '’Cold Pw ® Pcoid is then possible iff 

^a('’Cold ® Pwj fColdw) — 'fa(PCold ® Pwj '’Coldw) ^ ^ 0, (B7) 

where is the thermal state of the system at temperature Thoi of the surrounding bath. The quantity Ea{p, a) for a > 0 
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corresponds to a family of free energies defined in [5], which can be written in the form 

Fa{p,T) = ^[DaipWr)-\nZh], (B8) 

Ph 

where Z1 q,(p||t) are known as a-Renyi divergences. Sometimes we will use the short hand Fao ■= limo,_>.oo Fa- On occasion, 
we will refer to a particular transition as being possible/impossible according to the Fa free energy constraint. By this, we mean 
that for that particular value of a and transition, Eq. (B7) is satisfied/not satisfied. The a-Renyi divergences can be defined for 
arbitrary quantum states, giving us necessary (but insufficient) second laws for state transitions [5, 10]. However, since we are 
analyzing states which are diagonal in the same eigenbasis, the Renyi divergences can be simplified to 

Dc.{p\\t) = ^ (B9) 

i 

where pi^ qi are the eigenvalues of p and the state r. The cases a = 0 and a 1 are defined by continuity, namely 

Dq{p\\t) = lim 0„(p||r) =-In qi, 

Di{p\\t)= lim H„(p||t) = , 

I 

and we also define Doo as 

-Doo(p||r) = lim Oa(p||T) = Inmax —. 

a—>^00+ i Qi 

The quantity Di{p\\t) coincides with Z1 (p||t), as we have defined in Eq. (B3), and evaluated in Eq. (B4) for diagonal states. 
We will often use this convention. Eurthermore, since we are considering initial states which are block-diagonal in the energy 
eigenbasis, these generalized second laws are both necessary and sufficient conditions for state transformations. Therefore, in 
Section E2 a we will solve Eq. (B7) explicitly to find an expression for IText with the ultimate goal of evaluating the maximum 
efficiency in this regime. 

The reader should note that for both Section B 1 and B 2, the conditions for state transformation place upper bounds on the 
quantity H4xt- In particular, this allows us to express the maximum values Wext can take (such that the joint state transformation 
of cold bath and battery is possible) in terms of quantities related to the cold bath, and the error probability e. It is also worth 
comparing the conditions for state transformation in Section B 1 and B 2, which are stated in Eqs. (Bl) and (B7). In particular, 
Eq. (BI) is but a particular instance of Eq. (B7), and therefore the nanoscopic second laws always place a stronger upper bound 
on Wext compared to the macroscopic second law. 


(BIO) 

(BID 

(B12) 


C. Efficiency, maximum efficiency and how to evaluate it 

The central quantity of interest in this letter is the efficiency of heat engines. Since we have already introduced the notion 
of a heat engine in Section A, and the rules which govern the possibility of thermodynamical transitions of one cycle of a heat 
engine in Section B, it is timely to define the efficiency. After defining this quantity, we demonstrate how go about calculating its 
maximum value under different conditions, such as for perfect work, near perfect work, in both the macroscopic and nanoscopic 
regimes. This will prepare the scene for Sections D and E, where we evaluate the maximum efficiency more explicitly. 


1. Definition of efficiency and maximum efficiency 


As stated in the main text, the efficiency of a particular heat engine (recall that a heat engine is defined by its initial and final 
states PcoidHotMW: P&idHotMW as described in Section A) is defined as 


77 := 


M^ext 

AH ’ 


(Cl) 


where Wext is the amount of work extracted which is defined in Eq. (A7), and AH is the amount of mean energy drawn from 
the hot bath, namely AH := tr(77HotPHot) ~ tr(^HotPHot)’ where is the reduced state of the hot bath. 

Now, consider the set of conditions on state transformations given by Eq. (B7) for nanoscale systems. As discussed in Section 
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B, these conditions place a restriction on the range of values Wext can take. Therefore, for any fixed Pcoid’ define t7"™°(pcoid) 
as the maximum achievable efficiency as a function of the final state of the cold bath. More precisely, 

^““(Pcoid) (C2) 

= sup 7?(pcoid) subject to (Pw ® Pcm, Tqmw) Va > 0. (C3) 

Wex, 

In Eq. (C2), we have written the quantity in Eq. (Cl) as p = p(pcoid) remind ourselves of its explicit final cold bath state 
dependency. Therefore, the maximum efficiency will correspond to maximizing over the final state of the cold bath: 

Pmax = sup p”“°(pCold). (C4) 

t'cold^'S 

where S is the space of all quantum states in Tfcoid- By analyzing this quantity in Section E, we show that perfect work cannot 
be extracted. Therefore, when we calculate the maximization in Eq. (C4) we will consider near perfect work (see Def. 2). 


In the macro regime, we have to satisfy a less stringent requirement, namely the macroscopic second law of thermodynamics. 
And hence we have that for fixed p^oid’ P™^'^(Pcoid) maximum efficiency as a function of 


(Pcoid) = sup p(Pcoid) subject to 

Wex, 


-^(PColdMw) — ^(PColdMw) 


tr ( Ht PColdHotMW ) — tr ( iTt PColdHotMW ) ! 

where iJcoidHotMW is defined in Eq. (AI). Similarly to the nanoscale setting, the maximum efficiency is 

^max 


= sup r7”“(pcoid)- 


Pcold^'S 


(C5) 

(C6) 

(C7) 


We can also define the maximum quasi-static efficiencies for the macro and nano scale. The maximum efficiency of a quasi-static 
heat engine (see Def. 3), is 


= lim p"“°(r( 5 )), 

9^0+ 

(C8) 

= lim p”-(r( 5 )), 

3^0+ 

(C9) 


for the nanoscopic and macroscopic cases respectively. T{g) G "Hcoid is the thermal state with Hamiltonian ITcoid at temperature 
Pf = Pc — g and are defined in Eqs. (C2) and (C5) respectively. Since we can extract perfect and near perfect work 

in the macroscopic setting, we will derive the efficiency for both cases in Section D. 


2. Finding a simplified expression for the efficiency 


We can find a more useful expression for AH appearing in Eq. (Cl). This can be obtained by observing that since only energy 
preserving operations are allowed, we have 


tr( PColdHotMW ) — tr (-^t PColdHotMW ) I (CIO) 

where Ht = TTnot + TTcoid + Hyi -f T/w- Since the Hamiltonian does not contain interaction terms between these systems, the 

mean energy depends only on the reduced states of each system. Mathematically, it means that Eq. (CIO) can be written as 

tr(iTHotpSot) + h-(iTcoidPcoid) + ^AHupIi) + tr{Hwp%) = (Cl 1) 

tr(iTHotPHot) + tr(fTcoidPcoid) + ^(HuPm) + tr(iTwPw)- (C12) 


Also, note that since therefore tr(i7^MPM) = 'r(iTMPM)- This implies that we have 


AH = AC + AIE, 


(C13) 
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where 


AC := tr 


-ffcoldPcold 


- tr 


-ffColdTcold 


(C14) 


and 


AW := tr(fi^wPw) “ tr(iJwPw)- 


(C15) 


are the change in average energy of the cold bath and battery. We can thus write Eq. (Cl) as 




V = 


AW + AC 


(C16) 


Furthermore, from Eqs. (A5), (A6), (A7) and (C15), we have AW = (1 — s)Wexi, and hence we can write the inverse efficiency 
as 


V ^(Pcoid) = l-e + 


^^(PCold) 

4Eext(Pcold) 


(C17) 


where we have made explicit the dependency. We already know from the setting that is thermal. If is also a 
thermal state at some temperature /3 according to the cold bath Hamiltonian 7/coid, we will sometimes use the shorthand notation 
VW) for77(pcoid) and AW{/3), AC(/3) for AIE(pJ„iJ, AC(p^„,J respectively. 

In Section D, we will derive an expression for Wext and solve Eqs. (C5), (C7). In Section E, we will derive a new expression 
for Wext in the nanoscopic regime, and solve Eqs. (C2), (C4). 


D. Efficiency of a heat engine according to macroscopic thermodynamics 

In this section, we study the efficiency of the setup detailed in Section A under the constraints of macroscopic thermodynamics, 
as described in Section B 1. This implies that the Helmholtz free energy solely dictates whether Pcoidw Pcoidw i® possible. 
We find that in both cases of extracting perfect and near perfect work, 

(1) The maximum achievable efficiency is the Carnot efficiency. 

(2) The Carnot efficiency can be achieved for any cold bath Hamiltonian. 

(3) The Carnot efficiency is only achieved when the final state of the cold bath is thermal (according to a different temperature 

Tf)- 

(4) The Carnot efficiency is only achieved for quasi-static heat engines, meaning in the limit where 7/ > Tcoid- A technical 
definition of quasi-static heat engines can be found in Def. 3. Roughly speaking, this means that there is only infinitesimal 
change in the final temperature of the cold bath, compared to its original state. 

This section can be summarized as follows: in Section D 1, we first apply the macroscopic law of thermodynamics, namely 
the fact that Helmholtz free energy is non-increasing, to our heat engine setup. By making use of energy conservation, we can 
derive the amount of maximum extractable work as shown in Eq. (D4). Next, in Section D 2 we show that when considering 
the extraction of perfect work, we show the points (l)-(4) as stated above. In Section D 3, we show that points (l)-(4) hold also 
when considering near perfect work. 

The main results can be found in Theorem 1 and Lemma 6. One may think points (l)-(4) are obvious since it has long been 
known that the optimal achievable efficiency of a heat engine operating between two thermal baths is the Carnot efficiency, and 
that this efficiency can only be achieved quasi-statically. The motivations for proving these results here are two-fold. Firstly, 
this is a rigorous and mathematical proof of optimality, while usually one encounters arguments such as reversibility, or that 
the heat engine must remain in thermal equilibrium at all times during the working of the heat engine. Secondly, we will find 
later on at the nano/quantum scale that the Carnot efficiency can be achieved but observation (2) does not hold anymore. For 
these reasons, it is worthwhile proving that one can actually achieve points (l)-(4) in this setting for any cold bath Hamiltonian 
according to macroscopic thermodynamics. From a practical point of view, many of the technical results proved here will be 
needed in the proofs of Section E, where we derive results involving a more refined set of generalized free energies, which 
describes thermodynamic transitions for nanoscale quantum systems. 
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1. Maximum extractable work according to macroscopic law of thermodynamics 

Our first task is to find an expression for Wext in the macro regime. We do so by solving Eq. (Bl) for Wext such that 

— ^>S'(PcoldMw) — — ^>5'(PcoldMw)- (Dl) 

The entropy is an additive quantity under tensor product, meaning that S{pi ® P 2 ) = S{pi) + S{p 2 ) for any states pi,p 2 - 
Furthermore, since the joint Hamiltonian does not contain interaction terms, therefore the mean energy also depends only on the 
reduced states. In summary, both S and (H) are additive under a tensor product structure of Pcoi^mw PcoWmw described 
in Eqs. (A2) and (A3). This means one can rewrite Eq. (Dl) by expanding its terms, 

(-^Coid)pi_^j^ + + ('^w)pi^ - ^ ['S'(pLid) + 5 '(Pm) + 'S'(Pw)] ^ (D 2 ) 

(ffcold)pO^jj + (-^M)pO + ('^w)pO, - ^ ['S'(Pcold) + 'S'(Pm) + '5'(Pw)] I 

Furthermore, note that = p^, and therefore S'(p^), {Hm)p^ are common terms on both sides of Eq. (D2) which can be 
cancelled out. Furthermore, by our construction of the battery in Eqs. (A4)-(A7), we have that S{p^) = 0, <S'(p^) = AS' = 
h 2 (e) and (iTw)pO = and (H^w)p^ = E'^. Thus, Eq. (D2) can be simplified to 

W4xt + (-ffcold)pi^|_j — ^5'(pcold) < (^Cold)pC!^|^ — ^>S'(Pcold) + ^h 2 (e), (D3) 

where Wext has been defined in Eq. (A7). In other words, Wext < ^(Pcoid) ~ -^(Pcoid) + /) 7 ^ 2 (e)- 

We can also express Wext with the relative entropy instead, by using Eq. (B 6 ). We can apply this identity to Eq. (D3) whenever 
the initial and final states are diagonal in the energy eigenbasis. Note that the initial Pcoid a thermal state (of some temperature), 
and therefore diagonal in the energy eigenbasis. Since we start with a state ® p^ which is diagonal w.r.t. the Hamiltonian, 
and since catalytic thermal operations can never increase coherences between energy eigenstates (or in the macro setting, since 
we only demand mean energy conservation), we know that the final state Pcoid ® Pw diagonal in the energy eigenbasis. 

Therefore, Eq. (D3) can be rewritten w.r.t. the relative entropies as follows 

hkext < ^"(PCold) ~ ^(Pcold) + ^h 2 (e) = — [^'(Pcoldlkcold) ~ ^(PColdlkcold) +h 2 (e)] • (D4) 


2. Maximum efficiency for perfect work is Carnot efficiency 

In this section, we want to find the maximum efficiency according to Eqs. (Cl), (C5) and {Cl), for the case of e = 0 which 
implies h 2 (e) = 0. We do this by the following steps: 

1. Evalaute Wext- According to Eq. (D4), we know that 

Wext = -F(pcoid) — -F(Pcold) = ^ [-O(PColdlkcold) ~ -^(PColdIkcold)] i (D5) 

where recall that we have defined previously as the thermal state of system C with temperature Thoi- Note that here 
equality can be achieved because in macroscopic thermodynamics, satisfying the free energy constraint is a necessary and 
sufficient condition for the possibility of a state transformation. Note that since by construction the initial and final states 
of the battery are pure energy eigenstates, namely e = 0 and therefore 

Wext = AW. (D6) 

2. Write inverse maximum efficiency as optimization problem. By substituting the simplified expression for efficiency derived 
in Eq. (Cl 7) into Eq. (C7), we have 


7 


-1 

max 


= inf (77“’"'=)-^ 

fcold 


= 1 + inf 

T’colil 


AC 


(D7) 
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3. Maximize Wgxt given a fixed value of AC. This is done in Lemma 1, where we show that given a fixed AC, the final cold 
bath state that maximizes Wext is uniquely a thermal state, corresponding to a certain inverse temperature /3'. 

4. Show that 3) implies that efficiency is maximized by a thermal state of the cold bath. This is proven in Lemma 2. Therefore, 
this implies one only needs to optimize Eq. (D7) over one variable, i.e. fif, the final temperature of the cold bath. 

5. Show that the efficiency is strictly increasing with fif. This is done first by proving several identities, which are summarized 
in Corollary 1. Using these identities, we prove in Lemma 4 that the first derivative of efficiency w.r.t. /?/ is always positive 
over the range where Wext > 0. This leads us to conclude, in Theorem 1, that maximum efficiency is achieved in the limit 

> /3c, and evaluating the efficiency at this limit gives us the Carnot efficiency. 

Firstly, let us develop a technical Lemma 1, which concerns the unique solution towards maximizing Wext for a fixed AC. By 
applying Lemma 1, we show in Lemma 2 that the maximal efficiency is achieved when Pcoid ^ thermal state. The reader can 
easily find similar proofs in [29]. 

Lemma 1. Given any Hamiltonian HcoU, o corresponding thermal state of some temperature fih, eind a fixed initial state 

Pcoid’ consider the maximization over final states Pcey, 


max Wext (D8) 

t'cold 


over all states Pcoid which are diagonal in the energy eigenbasis, subject to the constraint that AC is a constant. Then the 
solution for is unique, and p is a thermal state according to the Hamiltonian TTcoid ot a certain temperature fi'. 


Proof. Firstly, from Eq. (C14) we see that the constraint AC being a constant, is the same as tr iLcoidPcoid being a constant. 
This is because they differ only by a constant term. On the other hand, from Eq. (Cl5) and (D6), we can see that Eq. (D8) is 
equal to 


rnax Wext = yr 

^Cold 


f^(PColdlkcold) 


minL>(pcoidll "^Cold) 
Pcq\6 


(D9) 


Since p^oid ^ both diagonal in the energy eigenbasis (Pcoid by the statement in the lemma, and r by it being a thermal 
state), one can evaluate the relative entropy by using Eq. (B3). Denote the eigenvalues of our variable p^oid be {pi}i, and the 
eigenvalues of the thermal state r to be {qi}i. We can then write the optimization problem as 


min 

{pi} 


'^Pi{\npi - Ing,); 
i 


where 



subject to ''^^piEi = c constant, and ''^^pi = 1. 

i i 


We can now employ techniques of Lagrange multipliers to solve this optimization. The constrained Lagrange equation is 


dL 

dpi 

dL 

d\ 

dL 

dp 


'^pfilnpi - Inq^) + X 

i 




(InPi - In gi + 1 + XEi + p) = 0, 


EiPi - c = 0. 

i 

-1 = 0 . 
i 


We find that the normalized solution is 


Pr = 



Zp, = 


(DIO) 

(Dll) 

(D12) 

(D13) 

(D14) 


and Pi are probabilities corresponding to the Boltzmann distribution, according to inverse temperature /3' = /3 + A. Depending 
on the mean energy constraint c and normalization condition, one can solve for the Lagrange multipliers A and p. With this we 
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conclude that the state p which maximizes £^(pcoidll''’) ^ thermal state, where its temperature is such that the constraint on 
mean energy is satisfied. □ 

Lemma 2. Consider the work extraction process described by the state transformation PcoidMW PcoidMW> tv/zere p^oid’ Pw 
and p\i have been described in Section A. Denote 'Hco\d os the Hilbert space of the cold bath. Then the maximal efficiency in 
Eq. (D7) is obtained for a final state of the cold bath Pcoid’ which is thermal: 


V 


-1 

max 


1 + inf 


AC 


(D15) 


where Sr the set of all thermal states (for /fcoid with any temperature T > 0) in "Hcoid- Furthermore, all non-thennal states do 
not achieve the maximum efficiency, i.e. 


7max 


< 1 + 


AC 


We, 


PcoM 


for any eS\Sr 


(D16) 


where S is the space of all quantum states in "Hcoid 

Proof First of all, note that without loss of generality we can always consider only diagonal states, as explained in the para¬ 
graph before Eq. (D4) that catalytic thermal operations do not increase coherences between energy eigenstates. We begin by 
substituting Eqs. (C14) and (D5) into Eq. (D7), and finding 


= 1-1- inf 

Z^Cold 


AC 


= 1 -f inf 


PhAC 


Pcoid ^l('^Coldll'^Cold) -^l(f^Coldll'^Cold) 


— 1 + /3?i 


^l('^Coldll'^Cold) -^l(f^Coldll'^Cold) 

olj.^3 ..N .A. 

P^oid tr(iFcoidPcoid) - tr(^^Coid'r^oid) 


-1 


(D17) 

(D18) 

(D19) 


In the last line of Eq. (D19), we see that only two terms depend on the maximization variable Pcoid- This means we can perform 
the maximization in two steps: 


1 ('^Cold 11 '^Cold ) 1 ( Pgp) j 11 Tcoid ) _ Bi^A) 

PcaXd tr(iFcoidPcoid) - tr(i7coidr^oid) ^>o ^ 


(D20) 


where B{A) is the optimal value of a separate minimization problem: 

B{A)= inf i9i(pi„idlkc^id) (D21) 

PcoldSf 

tr{ Cold Pcold ) ~ 'T ( Hcq\A ^Cold ) ~ ^ 

Erom Lemma 1, we know that the solution of the sub-minimization problem in Eq. (D21) has a unique form, namely Pcoid = '^Coid 
is a thermal state of some temperature /3j. Therefore, Eq. (D20) can be simplified to 

g^p -Pl('^Coldlkcold) ~ T>l(Pcoldlkcold) _ g^p -^>1 (Tgpjd|| — -Pi (fcoid Ikcold) (D22) 

pJoid ti-(^CoidPcoid) - tr(^Coidr^oid) pf tr(iJcoidT-^oid) “ tr(^CoidT^oid) 

Whats more, for every constant A, the function 


fix) 


1 ^ j- R, 

T^l('^Coldb'Cold) X 

\ ^ ^ Ph 

A 


(D23) 


is bijective in a: S R and thus due to the uniqueness of the sub-minimization problem in Eq. (D21), we conclude that for all 
non-thermal states Pcoid’ corresponding efficiency will be strictly less than that of Eq. (D19). Thus from Eq. (D22) and (D19) 
we conclude the lemma. □ 

After establishing Lemma 2, we can continue to solve the optimization problem in Eq. (D7) by only looking at final states 
which are thermal (according to some final temperature (3f which we optimize over). In the next Lemma 3 and Corollary 1, 
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we derive some useful and interesting identities. These identities concern quantities such as the derivatives of mean energy 
and entropy of the thermal state (with respect to inverse temperature), and relates them to the variance of energy. We later use 
them in Lemma 4 to prove that the Carnot efficiency can only be achieved for quasi-static heat engines. The reader can find 
similar proofs in any standard thermodynamic textbook (For example in Sections 6.5, 6. of [23]), but we derive them here for 
completeness. 

Lemma 3. For any cold bath Hamiltonian fTcoid. consider the thermal state Tp = with inverse temperature /3. 

Define (Ffcoid)/? = tr(-ffcoidL3)’ mean energy and entropy ofrp. Then the following identities 

hold: 


d{Hco\d)p ,£r \ 

-= -var(Tfcoid)/3 

—— = -/3 . var(iTcoid)/3, 


(D24) 

(D25) 


where var(Tfcoid )/3 = {HcoXd^P “ (^^Coid)| is the variance of energy for rp. 

Proof. Intuitively we know that the expectation value of energy increases as temperature increases (or as the inverse temperature 
decreases). More precisely, consider the probabilities of Tp for each energy level of the Hamiltonian Ei, 


Pi = 


dpi 

dp 


3-/3Si 


where Zp = ( 


-PEi 


1 


-EiC 


-PE, 


Za- 


dZ, 


dp 


I . ^-PE, 


-PiE^ 


1 dZp 

^pliP 


Pi = -PiEi + Pi{Hco\d) p- 


(D26) 


The last equality holds because of the following identity: 


= (^Co.d)/3. 


(D27) 


Therefore, we have 


d{Hcoid)p _ ^ d{Hco\d)p dpi 


— ^ ^ Ei • PiEi Pi(ffcold)/3 


(D28) 


dp ^ dpi dp 

i L 

= -{El.xd)p + (^coid)^ = -var(i7coid)/3. (D29) 

On the other hand, similarly, one can prove the second identity by writing down the expression of entropy for the thermal state. 


S(/3) = - E 


J-SE. J-SE. 

• In —— 


e-PEi ^ ^ 






2 « 


= P{Hco\d)p + ln2^. 


(D30) 


Therefore, the derivative of S{P) w.r.t. P is 


^ ^ (He* + ^ ^ ^ ^ ■ var(Ac„)E. 


□ 


By using Lemma 3 in a special case, we obtain the following corollary: 
Corollary 1. Given any Hamiltonian Tfcoid. consider the quantities 


AC{Pf) = tr{HcoidEpf) — t.r{HcoidEpPj = {Hco\d)Pj — {Hcoid)p„ 


(D32) 


and 


Wext{/3/) = F{TpJ - E{Tpj) = ^ [D{TpP\TpP) - D(Tpj\\TpJ\ , 


(D33) 
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where Tp corresponds to the thermal state defined by i/coid at inverse temperature (3. Then 

dACiPf) , 

= -var(i?coid)/3, 

dWsniPf) Ph-Pf ^ 

Proof. For AC'(/3/), it is straightforward from Lemma 3 that 

dAC{/3f) d{Hcou)i3} 

-^PT “ “ -v.r(ffco„)s,. 

On the other hand, AfF(/3/) can be simplified by substituting Eq. (B6) into Eq. (D33), 

Wext(/3/) = Fijfif) - F{Tfij) = {Hco\a)p^ - {Hco\A)fij - ^ [5'(t;3J - . 

With this, we can evaluate the derivative 


dWUPf) 

dl3f 


d{Hcoid)0f 1 dSiTfjj) 


dPf 


Ph dPf 


= var(iFcoid)/3y - ^var(iJcoid)/3 


^-Pf 

Ph 


var(iFcoid)/3j. ■ 


(D34) 

(D35) 


(D36) 


(D37) 


The second equality is obtained by Lemma 3 for 


d{Hca\d) 0 . 
Wf 


and the third by grouping common factors together. 


□ 


In the next step, by using Corollary 1, we show that the optimal efficiency is achieved only in the quasi-static limit, i.e. in the 
limit Pf —>• Pc- 

Lemma 4. Evaluate the efficiency expressed in Eq. (Cl 7) /or the situation where the final state of the cold bath is a thermal 
state at inverse temperature Py: 


viPf) 


WM) 

AciPf)+wMy 


(D38) 


Then for all Pf < Pc, ^ > 0. 


Proof. To prove this, we show that 


df£f_ 

dPf 


< 0, where rj ^ 


1 -f ^r-. Evaluating the derivative of ry ^ w.r.t. Pf, we obtain 


dr] ^ 

w 




dACiPf) dWcAPf),^ 

-^pt . 

-W/xt - AC 

Ph 


var(i7coid)/3/ 
var(77^Coid)/3/ 
var(i7coid)/3/ Pf 


AC + ^ 




Ph 


AC - -lS(Tp,) - S(TpJ] 


(D39) 

(D40) 

(D41) 

(D42) 

(D43) 


The first equality is obtained by invoking the chain rule of differentiation. The second equality is obtained by substituting 
evaluated earlier in Corollary 1 . The third equality is obtained by expressing Wext according to Eq. (D37), plus 
recognizing that {Hcoitpr/B^ — {Hcoid)Tfi^ = AC. The last inequality is obtained, simply by taking out a common term pf /Ph. 
































We then make the following observations: 

1) The factor 
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PhWl, 


> 0 , 


2) The variance of energy for any positive temperature 


var(7Tcoid)/3j > 0, 


(D44) 


(D45) 


3) and the last term AC — — S{tis^)] can be written as F{Tpj) — F{Tp^), where F is the free energy of a system w.r.t. 

a bath with inverse temperature f3f. But then, since is the thermal state with the same inverse temperature, this means that 
Tpj is the unique state that minimizes free energy. Therefore, F{Tp^) — F{Tp^) > 0 for any 

□ 


From Lemma 2 and Lemma 4, we conclude that the maximization of efficiency for any Hamiltonian F[ happens for a final 
state which is thermal, and the greater its inverse temperature /3y, the higher efficiency is. With these lemmas we can now prove 
the main result of this section (Theorem 1). 

In the next theorem, we evaluate the efficiency at the limit /3/ —/3“, and show that it corresponds to the Carnot efficiency. 


Theorem 1 (Carnot Efficiency). Consider all heat engines which extract perfect work (see Definition 1 ). Then according to the 
macroscopic second law of thermodynamics, the maximum achievable efficiency (see Eq. (Cl)) is the Carnot efficiency 


= 1 - 


I3h 

Pc' 


(D46) 


It can be obtained for all cold bath Hamiltonians Ffcoid. but only/or quasi-static heat engines (as defined in Def. 3 and Eq. (C9) 
for quasi-static maximum efficiency), where an infinitesimal amount of work is extracted. 


Proof. From Eq. (C7), we have an expression for the optimal efficiency in terms of a maximization over final cold bath states 
Pcoid € S. By Lemma 2, we know that the optimal solution is obtained only for thermal states. Subsequently, by Lemma 4, it 
is shown that when the final cold bath is of temperature (3f, the corresponding efficiency is strictly increasing w.r.t. jdf. Also 
note that since by definition Wext > 0, this implies that /3/ < /3c- Intuitively, this is because heat cannot flow from a cold to hot 
system without any work input. One can also see this mathematically, by showing that for any P > /3h, 

£ - var(iTcoid)/3 > 0. (D47) 


dFjTp) 

dp 


pL 

dp 


(-ffcold)/3 - 


This implies that if Pf > Pc > Ph, then F(Pf ) > F(Pc), and according to Eq. (D33) Wext A 0. Therefore, the optimal 
efficiency must be achieved only when the Pcoid ts a thermal state whose inverse temperature /3/ approaches Pc from below. Let 
Pf = Pc — g, where g > 0. Then Thus we have 


p-i = \un ir,--)-\Pc - g), - fl) = 1 + 

9^0+ 


AC 

w;^ 


Pcold~'’'{0c-9) 


(D48) 


Since as p —0+, both the numerator and denominator vanish, we can evaluate this limit by first applying L’Hospital rule, the 
chain rule for derivatives (for any function F, ^ then Corollary 1 to obtain 


lim 

g->o+ 


AC 


lim 

g^0+ 


dAC 

dg 

'dW^ 

dg 


lim 


dAC 

ddf 

dWcn 

ddf 


Ph 

Pc- Ph 


This implies that 


7 


-1 

max 


lim (77”-)-i{^c 

3^0+ 


7) = 1 + 


Ph 

Pc - Ph 


Pc 

Pc - Ph 


(D49) 


and hence r/max = 1 - f;- 


□ 
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3. Maximum efficiency for near perfect work is still Carnot efficiency 


In this section, we show that even while allowing a non-zero failure probability e > 0 in the near perfect work scenario, the 
maximum achievable efficiency is still the Carnot efficiency. It is worth noting that this result is also important later, as an upper 
bound to maximum efficiency in the nanoscopic regime. We first prove it in Lemma 5 for the case where the final state of the 
battery is fixed as in Eq. (A6). Then later, we show in Lemma 6 that Carnot is still the maximum, even if we allow a more 
general final battery state. Before we present the proof, it is useful for the reader to recall the definition of near perfect work 
(Def. 2) and quasi-static heat engines (Def. 3). 

Lemma 5. Consider all heat engines which extract near perfect work (see Def. 2). Then according to the macroscopic second 
law of thermodyanmics, the maximum efficiency of a heat engine, p^nax is the Carnot efficiency 

Vmax = sup ?7™“(Pc) = ^ (^50) 

P^eS Pc 

and the supremum can only be achieved for quasi-static heat engines (see Def. (3) and Eq. (C9)j. 

Proof. The ideas in this proof are very similar to that of Section D 2, and the main complication comes from proving that even 
if we allow e > 0, as long as AS'/lLext is arbitrarily small, the maximum efficiency cannot surpass the Carnot efficiency. 

Let us begin by establishing the relevant quantities for near perfect work extraction. The amount of work extractable from the 
heat engine, when we have a probability of failure, according to the standard free energy can be obtained by solving Eq. (D5). 
We thus have that the maximum Wext is 

W^ext = pf\l - £)-' [D{tp^\tp,) - + A5] , (D51) 


where AS” is defined in Eq. (A9). 

Before we continue with the analysis, we will note a trivial consequence of Eq. (D51). Condition 1) in Def 2 implies that 
(1 — e)“^ is upper bounded. The terms in square brackets in Eq. (D51) are also clearly upper bounded for finite Pc^Ph- Hence 
Wext is bounded from above. AS' is solely a function of e and only approaches zero in the limits e O"*", e —1“; and e —1“ 
is forbidden by 1) in Def 2. Thus if 1) and 2) in Def 2 are satished. 


lim ^ = 0. 
£—>■0+ Wext 


(D52) 


In turn, if Eq. (D52) is satished, then we have near perfect work by Def. 2. Thus Eq. (D52) is satished iff we have near perfect 
work. We will use this result later in the proof. 

Extracting a positive amount of near perfect work implies that we can rule out all states p}. such that D{Tp^ \ \ Tp^ ) < D{pq\\ rp,, ) 
from the analysis. This can be proven by contradiction; if D(TpJ\Tp^) < i9(pp||T/3^), thenfromEq. (D51) /3?tWext < AS/(1—e) 
and together with 2) in Def 2 this would imply 

AS 

0</3/i(l-e) < <p. (D53) 

^^ext 


However, since from 1) Def. 2 we have e < I, Eq. (D53) cannot be satished for all p > 0, leading to a contradiction. 

Erom Eq. (D7) we have 

AC 

7max = 1 - £ + inf = (1 - e) • 

Wext 

where AC = AC(pq) and is dehned in Eq. (C14). 

Eirstly, let us show that with a similar analysis as shown in Lemma 2, the maximum efficiency occurs when p^ is a thermal 
state. Erom Eq. (D54), we have 


1 + 


PhAC 


D{Tp\\Tp.) - lirflj -f AS 


(D54) 


V 


-1 

max 


(l-e) 


(1-e) 


1 I • r _ _ _ _ __ 

^ i9(r;3:iKr--“C(pi IIt;,;) + AsJ 

1 ^ ’ f ^ 

^^^XloDirpM-BiA) + AS_ 


(D55) 


(D56) 
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where 


B{A)= inf D{pI,\\tp,). (D57) 

PcG5 

t»’(ffcoldPc)-tt(-WcoldT/3„ ) = A 

We can split this minimization problem to Eqs. (D56) and (D57) because D{TpJ\\Tp^) and A5' do not depend on the variable 
Furthermore, when Pcoid ^ thermal state of inverse temperature /3/, we have seen in the beginning of the proof in Theorem 1 
that for IText > 0, /?/ < /3c- This implies that the variable A = AC = tr(iTcoid'rS/) ~ tr(37^coidT/3c) > 0- 

By Lemma 1, for any fixed A > 0 we conclude that the infimum in Eq. (D57) is achieved uniquely when is a thermal state. 
Therefore, our optimization problem is simplified to optimization over final temperatures /3/ (or p = /3c — /3/), 

7mL = (!-£)• 


AC 

AOO 


AS 


(D58) 


Consider cases of /3/, where D{TpJ\\Tp^) — I3(Tpj. ||rp^) is non-vanishing (finite), which are non quasi-static. Note that this 
always corresponds to extracting near perfect work, since when e -A O’*", we have e, AS -A 0 and these contributions dissapear 
from Eq. (D58). However, by Lemma 2 we also know that the infimum over /3/ occurs uniquely at the quasi-static limit, when 
p —O'*'. This means that for all non quasi-static cases, Carnot efficiency cannot be achieved. 


What remains, is then to consider the quasi-static heat engine, namely the limit p —O'*'. Extracting near perfect work in this 
case corresponds to requiring that limg_^Q+ = 0, where e = e{g) and limg_^o+ ^{ 9 ) = 0- Equivalently 


lim 

g->-0+ 




AS 


= 00 . 


(D59) 


Substituting Eq. (D51) into Eq. (D59), 


lim (1 - s{g)) ^ 1 

g->- 0 + 


-P(T'/3clkpJ --P(rp^||rpJ 
AS- ^ 


= 00 


which implies that lim 

g-*'0+ 


D{TpJ\TfiJ - -D(rpJ|rpJ 
AS 


00 , or equivalently. 


(D60) 


lim lim 


AS- 


e^o+ g^o+ -D(rpJ|rpJ -i9(rp^||rpj 


= 0 . 


(D61) 


Finally, we evaluate the inverse efficiency at the quasi-static limit, 

AC 


V = lirn (l-£(5')) 

g^0+ 

= l-\- Ph lim 

g-^■o^ 

= l-\- Ph lim 

g-i-oH 

= l-\- Ph lim 
= 1 - 


1 + /3/t 


- DiTpfWrpJ + AS _ 
AC 


9^0+ £i(TpJlTpJ-i9(Tp^||rpJ-f AS* 
AC 


g^0+ [-D(TpJ|rpJ-D( tpj.||tpJ] 
dAC{Tpj)/dg 
gAo+ di9(Tpj.||rpJ/dp 

Ph 


1 + 


AS* 


D{TpP\Tlj^) - D{T0j\\TfiP) 


-1 


Ph-Pc 


(D62) 

(D63) 

(D64) 

(D65) 

(D66) 


where from Eq. (D64) to (D65), we make use of Eq. (D61) : the second term within the limit is simply 1, and the first term 
depends only on p, which we can obtain Eq. (D65) by invoking the L’Hospital rule. The last equality in Eq. (D66) follows 
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directly from the identities we derived for and in Corollary 1 , 


dAC 

dg 

dPjTpfWTfiJ 

dg 


^ -var(i?coid)/3, 
dPjTp^WTf,^) ^ dWext 

d/3/ dp f 


iPh - /3/)var(fdcoid)/3y, 


while in the limit g —>■ 0, Pf = Pc- 

Finally, we now see that the quasi-static efficiency is 


I' Ph- Pc- Ph \ ^ _ Pc - Ph _ ^ Ph 
\ Ph- Pc ) Pc Pc 


(D67) 

(D68) 


(D69) 


which is exactly the Carnot efficiency. 


□ 


Later, in Section F2 we will need Lemma 2 to hold in a more general scenario, i.e. instead of the final battery state being 
= (1 — e) \Ek){Ek Iw + s \ddj){Ej Iw’ want to allow the final battery state to be any energy block-diagonal state with trace 
distance e. Next we state and prove this generalized lemma. 


Lemma 6. Consider all heat engines which extract near perfect work (see Definition 2), but allowing for any final battery state 
with a trace distance e to the ideal final pure state \Ek){Ek 1^. Then according to the macroscopic second law of thermodynamics, 
the maximum efficiency of a heat engine, rjmax is the Carnot efficiency 


^max — 


sup r] 
p(.&s 


=(4) = i-^. 


(D70) 


and the supremum is only achieved for quasi-static heat engines (see Def (3) and Eq. (C9)). 

Proof Firstly, let us note that since the initial state Pcoidw is energy block-diagonal, the final state has to also 

be block-diagonal. Therefore, given the product structure between the cold bath and battery, it is sufficient to consider the case 
when the final battery state is energy block-diagonal. Next, let us note that any final state which is energy block-diagonal, 
and has trace distance e with \Ek){Ek\^ can be written as. 


Pw = (1 - e) \Ek){Ek\^ + ep’w'^, where p’^"'' = |Si)(F^i|w > ^Pi = 1 andp^ = 0. (D71) 

i i 

Next, one can calculate Wext given by the standard free energy condition, i.e. 

+ E(p%) > F(pi„iJ + F(p}^). (D72) 

Using the identity E{p) = tr(idp) — P~^S(p), we have that 

F(tpJ + E, > F(phcu) + (1 - e)Ek + £tr(idwp'w') - PP^S{pl). (D73) 

Substituting Wext = Ek — Ej, and rearranging terms, we have 

(1 - e)H^ex. < F(rpJ - F(pi„u) + Pp^AS - e[tr(7dwp';"") - E,]. (D74) 


Finally, by using the identity (in Eq. (B6)) that F{p) = Pf^ ^[-D(p||rp^) — InZp^], the maximum amount of extractable work is 
given by 


Wext = (1 - e)-^pp^ ■ [i7(rpjlrpj - d7(pLidlkpJ + (D75) 

where E = tr(7d^wpw"'^) — Ej. 

Following the steps in Lemma 5, in particular the derivations in Eq. (D55) and (D56), we have 


V 


-1 

max 


(i-s). 


1 + Ph inf 

AOO 


_AC_ 

■D(rpJ|TpJ - Li(rp^||TpJ -f AS'-ei3 


(D76) 
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To show Eq. (D76) gives the Carnot efficiency, we show that 1) for non quasi-static cases where /?/ < j3c, Carnot efficiency is 
not attained, and 2) in the quasi-static limit, Carnot efficiency is attained. 

Let us first consider the case of extracting a non-vanishing amount of near perfect work, i.e. for all cases where Pf < Pc- 
Then near perfect work, by Def. 2, corresponds to the limit e —> 0, 


p ^ = lim (1 — e) 

S —^0 


— 1 + /3ti 


1 + /3ti 


AC 


- 0 (r/j^||T/ 3 j + AS - eE 


AC 


D{TfiJ\T0j - D{Tpf\\Tfi^)' 


(D77) 

(D78) 


In this limit, all terms involving e vanish, and the inverse efficiency has the same expression as the efficiency for perfect work. 
We already know from Lemma 4 that the infimum over Pf cannot be obtained in this regime, since the inverse efficiency is 
strictly decreasing with Pf. 

Therefore, again we are left with analyzing the quasi-static limit for this problem. Lollowing the derivation in Eq. (D64) for 
the quasi-static limit, we obtain 


»?max [D{TfiP\Tfi^) - D{Tpj\\Tfiy\ 


( AS-eE \ 

\ D{TfiP\Ti3^) - D{Tfif\\Ti3^) j 


(D79) 


where e = e{g) and note that requiring near perfect work implies that 

AS 

lim ----- ^ ^ = 0. (D80) 

g^0+ D{Ti3j\TfiJ - D{T0f\\TfsJ 

Next, we observe the relationship between e and AS, in the regime where e is small. Given any e > 0 denoting the trace distance 
d{p^, \Ek){Ek\yfP} = e, the smallest amount of entropy that can be produced corresponds to AS = /12(e). This is because if 
we try to distribute the weight e over more energy eigenvalues, then by majorization the entropy only increases. But we also 
know that e < /12(e) for small values of e, in particular over the regime e € [0, |]. Therefore, we have that in this regime, 
e < /12(e) < AS holds. Therefore, we also know that 


lim 


eE 


= 0 , 


(D81) 


g^0+ D{TpP\Tfi^) - D{Tljj\\Tp^) 

where e = e{g). Plugging Eqns. (D80) and (D81) into Eq. (D79), we have that the quasi-static efficiency is p = 1 — 4^. □ 


E. Efficiency of a nanoscopic quantum heat engine 

In this section, we will be applying the conditions for state transitions for nanoscale systems, as detailed in Section B 2. The 
reader will see that due to these extra constraints from the generalized free energies, the fundamental limitations on efficiency 
will differ greatly from those observed in Section D. 

Eirstly, in Section E 1, we show that the extraction of a positive amount of perfect work is impossible using the setup. In 
Section E 2, we show that this can be resolved by considering near perfect work instead. Then we hnd that 

(1) The maximum achievable efficiency is still the Carnot efficiency. This is proven in Section E2b. 

(2) However, the Carnot efficiency cannot be achieved for all cold bath Hamiltonians. This is our main result, which is 
stated in Theorem 2, found in Section E 2 f. The results in Section E 2 d and E 2 e are more technical proofs, that pave 
the way for deriving this main result. 

(3) The Carnot efficiency is only achieved when the final state of the cold bath is thermal (according to a different temperature 
Tf). This is proven in Section E 2 b. 

(4) The Carnot efficiency is only achieved for quasi-static heat engines (see Def. 3), meaning in the limit where Tf ^ Tcoid- 
This is proven in Section E2b. 
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1. Impossibility of extracting perfect work 


We will first show that with the general setup as described in Section A, no perfect work can ever be extracted. By this we 
mean that whenever e as defined in Eq.(A6) equals zero, then for any value of Wext > 0, and for any final state Pcoid’ '^^e 
transition |i?j)(Sj|w<8> T^oid \Ek){Ek\w^ Pcoid possible. Intuitively speaking, this occurs because the cold bath is 

initially in a state of full rank. Since thermal operations cannot decrease the rank of the system, therefore the final state of the 
cold bath must also be of full rank. By directly solving Eq. (B7), we find that the amount of extractable work satisfies 

Wext < kTiiot inf [-Da(7'Coldlkcold) ~ (Pcold Ikcold)] ) (Ef) 

Oi>(J 


where is the thermal state of the cold bath (according to the cold bath Hamiltonian Hcoid)^ at temperature Thoi (since the 
surrounding hot bath is of temperature Thoi)- However from Eq. (B8), -Do(tcoidll'tcoid) = ^o(Pcoidll'tcoid)' Therefore according 
to Eq. (El), the amount of work extractable satisfies Wext < 0. 

We phrase this with more rigor in the following Lemmas 7 and 8, which proves that for perfect work, lEext > 0 is impossible. 
The proof holds for general initial states Pcoid of full rank, in particular, they need not even be diagonal in the energy eigenbasis. 

Lemma 7. For any Wext > 0. consider the Flamiltonian given by Eq. (A4). Then for any inverse temperature Ph > 0, the 
thermal state tw = - 1 a , satisfies 

tr(e“Pfi"w) •' 


tl' [{\Ej){Ej\^N — \Ek){Ek\w) Tw] > 0. 


(E2) 


Proof. Eollows directly from the definitions. Since Wext > 0, we know that E, 

>0. 


tr(e 


< E'^ . Evaluating the quantity above gives 

□ 


Lemma 8. Consider any general quantum state Pcoid ^iffull rank. Then for any Pcoid> ^ke transition from Pcoid® Pw Pcoid® Pw 

is not possible via catalytic thermal operations if 


tr 


(np° - 


Tw 


> 0 , 


(E3) 


where Hp is the projector onto the support of state p, and Tw is the thermal state of the working body at the initial hot bath 
temperature. 


Proof. One can show this by invoking the quantum second law for a = 0 [5], which says that if pin —pout is possible via 
catalytic thermal operations, then 


DoipinWr) > Do{pout\\T), (E4) 

where r is the thermal state of the system at bath temperature, and 

I?o(p||cr) = lim —- lntr[p“CT^““] = — IntrppCr], (E5) 

a—>- 0 + O. — 1 

is defined for arbitrary quantum states p, a. Applying this law with pin = p^ 0 Pcoid und pout = Pw ® Pcoid’ arrive at 

T^o(Pwll''w) ~ f^o(Pwlkw) — f^o(Pcoldlkcold) ~ T)o(Pcoldlkcold)i (E6) 

where Tq^^^ and are thermal states of the cold bath and battery at the temperature of surrounding hot bath (Thoi) respectively. 
Since p^oid huve full rank, and since is normalized, therefore according to Eq. (E5), -Do(Pcoidll'^Coid) ~ *-*■ Eurthermore, 
since the a—Renyi divergence Dq is non-negative, therefore the r.h.s. of Eq. (E6) is lower bounded by 0. Thus, we have 

tr[(np“ - np^)'^w] < 0. (E7) 

Since this is a necessary condition for state transformations, we arrive at the conclusion that; when Eq. (E7) is violated, state 
transformations are not possible. But from Lemma 7, any type of perfect work extraction violates Eq. (E7). Therefore, in this 
setting, perfect work extraction is always impossible. 

□ 

To summarize. Lemma 8 implies that if the initial state of the cold bath is thermal, and therefore of full rank, then any work 
extraction scheme via thermal operations bringing p^ = |j)(j|wtop^ = |A:) (fc|w where Wext = E'^ — > 0 is not possible. 
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In general, we see that if IIpO 7 ^ 11^^, then when transition to is possible, transition to p^ is not. Consequentially, 
we will have to consider near perfect work at the nano regime. 


2. Efficiency for extracting near perfect work 

As we have just seen in the previous Section E 1, we cannot extract perfect work. Due to the impossibility result, we consider 
the relaxation of extracting near perfect work in the nanoscale setting. 

• We begin by evaluating the expression for efficiency according to the nanoscopic laws of thermodynamics, given a final 
state of the cold bath, and comparing it to the expression according to macroscopic laws of thermodynamics. This is 
done in Sections E2a and E2b, and the relation between two efficiencies are summarized in Eq. (E14). Since the 
nanoscopic efficiency is always smaller than the macroscopic efficiency, which attains Carnot efficiency only in the quasi¬ 
static limit, it will be possible only to attain Carnot efficiency in the quasi-static limit, when considering nanoscopic laws 
of thermodynamics. 

• We analyze the quasi-static regime, focusing on the special case where the cold bath consists of n qubits. Since the quasi¬ 
static limit corresponds to the case of small p > 0 , and e also has to be arbitrarily small for near perfect work extraction, 
we perform Taylor expansion of the analytical expressions for and AC w.r.t. g and e. This is done in Section E2c. 

• In Section E2d, we identify how to choose e{g) such that it corresponds to drawing near perfect work in the quasi-static 
limit. We first begin by observing that any continuous function e{g) that vanishes in the limit p —> 0 can be characterize 
with a real-valued parameter R that determine how quickly e goes to zero. This is shown in Lemma 11. In Lemma 12, we 
show that near perfect work is drawn only if R G [0,1]. 

• Lemma 12 gives us the analytical expression and minimization range in order to evaluate Wext, according to Eq. (E 66 ). In 
Section E 2 e, we show how one can evaluate this optimization problem, by comparing the stationary points and endpoints 
of the function that gives the leading term in Eq. (E 66 ). Lemma 13 proves a technical property of the first derivative 
of this function. Using it, we prove in Lemma 14 that one can always choose e(p) with some R < 1 such that the infimum 
of is obtained at either a = R or a ^ 00 . 

• Einally, in Section E 2 f, we use the results in Section E 2 e regarding the evaluation of Wext to find the efficiency in the 
quasi-static limit. 


a. An explicit expression for Wext 

Our first task is to work out an explicit expression for Wext depending on the initial and final states of the cold bath, e and hot 
bath (inverse) temperature /3^. Such as expression is found by applying the generalized second laws as detailed in Section B 2. 

Lemma 9. Consider the transition 


T-Coid Pw ^ Pcoid ® Pw with e > 0. (E 8 ) 

where and p\^ are defined in Eqs. (A5), (A 6 ) respectively. Let Wext denote the maximum possible value such that Eq. (E 8 ) 
is possible via catalytic thermal operations, with a thermal bath of inverse temperature fih- Tef fic > fih- Then the final state 
pLid = Cold is block-diagonal in the energy eigenbasis, and 


Wext = inf W„, 

Q :>0 

W„ = —-[ln(A-e“)-aln(l-e)], 

Ph{a - 1 ) 


(E9) 

(ElO) 

(Ell) 


where pi = 




= 




and p' are the probability amplitudes of state Pcoid ’^hen written in the energy eigenbasis of 


Tfcoid- The quantities Wi and defined by taking the limit a —>■ 1, -foo respectively. 
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Proof. Eq. (B7) is necessary and sufficient for Eq. (E8) to be satisfied. We can apply the additivity property of the Renyi 
divergence, to Eq. (B7) to find 


Da{p%\\T^N) + Da{Tp^\\Tp^) > £'a(Pwlkw) + -Da (PMd 1^/3 J, (E12) 

where is the thermal state with Hamiltonian iJw at inverse temperature Ph- We define Wa to be the value of E'^ — Ej^ 
that satisfies Eq. (E12) with equality. A straightforward manipulation of these equations gives the expression for Wa- Then 
W4xt = infQ,>o Wq is the maximum value that satisfies the inequalities Eq. (E12) for all a > 0. □ 


As we will see later there exist Pcoid ®tich that, Wext given by Eq. (E9), has a solution (i.e. Wext > 0) for any s > 0. 

We can use this to write down an explicit solution to the maximization problem Eq. (C2). Using Eqs. (C2), (C17) and Lemma 
9, we conclude 


_nano /I \ 

V (Pcold) 


1 — £ 


^^(PCold) \ 

inf„>o J 


(E13) 


where Wa is given by Eqs. (ElO), (El 1) and recall AC can be found in Eq. (C14). Erom Eqs. (E13), (ElO), (El 1), we see that 
the optimization problem sup^i ?7"™°(Pcoid) i® ^^11 a formidable task. In the next section, see will show that we can use the 
results from Section D, to drastically simplify the problem. 


b. An upper bound for the efficiency 

Before moving on to solving the nanoscale efficiency explicitly, we will first use the results of Section D 3 to find upper 
bounds for the efficiency in the nanoscale regime, in the context of extracting near perfect work (Def. 2). 

Recall how we have discussed in comparing Sections B 1 and B 2, that the solution for the family of free entropies Fa, in the 
case of El is simply the Helmholts free energy. Therefore, from Lemma 9, it follows that Wi is simply the maximum amount of 
extractable work according to Eq. (Bl). Erom Eqs. (C5), (C17), 

0 -^" to ) 

One can now compare Eq. (E14) with Eq. (E13), and note that for any p^oid ^ have Wi{pQ^^f) > infa>o Wa{Pco\i)- 

Therefore, we conclude that for any Pcoid ^ 

^““(Pcoid) < 7”“^(Pcoid)- (E15) 

Eq. (E15) in conjunction with Lemma 5 has an important consequence. Namely, 

sup ??"“°(Pcold) 

Pcom^‘5 

< \ — Ph/Pc if the state Pcoid '^hat solves the supremum is that of a quasi-static heat engine, 

< 1 — Ph/Pc if the state Pcoid *^hat solves the supremum is not that of a quasi-static heat engine 

This tells us that if we cannot achieve the Carnot efficiency for a quasi-static heat engine, we can never achieve it, and can only 

achieve a lower efficiency. We therefore will only consider the quasi-static regime in the rest of Section E. 


(E16) 

(E17) 


c. Evaluating near perfect work in the quasi-static heat engine 

In light of the results from the previous section, we will now calculate the near perfect work Wext for quasi-static heat engines, 
i.e. the case where e, g <C 1. Specifically, we make the following assumption about the cold bath Hamiltonian: 

(A.5) The Hamiltonian is taken to be of n qubits: 

n 

HccU = ^ ® He,k ® where He,k = Ek\Ek){Ek\, 


(E18) 




and E'fc > 0 is the energy gap of the fc-th qubit. 

The tensor product structure in Assumption (A.5) allows us to simplify Pcoid’ 
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n 



2 = 1 


(E19) 


where is the thermal state of ith qubit Hamiltonian Hi ^ at inverse temperature /3c. For the simplicity of following proofs, 
we present them in the special case of identical qubits, i.e. that Ei = E for all 1 < i < n. This means Eq. (E19) can be reduced 
to 


0 _ 

PCold — ■ 


Furthermore, since we consider quasi-static heat engines, the output state is 


„1 _ T-®" 

PCold — > 


(E20) 


(E21) 


with Pf = Pc — g ,where 0 <5< 1. Eq. (El8) together with Eq. (E21) allows us to further simplify Eq. (Ell) to 


A = 


fZMi:!] 


(E22) 


where pi = 


Hamiltonian He- The proof follows along the same lines as the proof to Lemma 9, but now noting that in Eq. (E12) we 
can replace /^^(rcoidll'r/j;,) and Da{pcc\d\\'^Ph) "'i^h nDa{TpP\TpfP) and nDa{Tpj\\Tpf^) respectively. This follows from the 
additivity property of the Renyi divergences. After proving the special case of identical qubits, we show in Theorem 2 that it can 
be extended to non-identical qubits as generally described by Assumption (A.5). 




Qi = 


C-l^hEi 


are the probabilities of thermal states (different temperatures) for the qubit 


Since we are dealing with near perfect work and quasi-static heat engines, both p > 0 and £ > 0 are inhnitesimally small. 
Thus with the goal in find of finding a solution for Wext from Eqs. (E9), (ElO), and (E22); we will proceed to find an expansion 
of Wa for small e and g. 


i) The expansion of A in a quasi-static heat engine 


To simplify our calculations of Wext, especially that of efficiency, it is important to express A in Eq. (E22) in terms of its first 
order expansion w.r.t. the parameter g. Recall that this parameter g = Pc — Pf is the difference of inverse temperature between 
the initial and hnal state of the cold bath. 

Firstly, note that for any integer n, the expression in Eq. (E22) evaluates to A|g=o = 1- This is because at g = 0, Pf = Pc and 
therefore the probabilities are identical. To obtain an approximation in the regime 0 < p <C 1, we derive 


dA 

dg 






dg 


= —an A 



Y^pTql-’^iE, 




(E23) 

(E24) 


The hrst inequality holds by noticing that only the probabilities p[ depend on g, which means only the denominator in Eq. (E22) 
is differentiated, using the chain rule 


dA{{p'i\) 

dg 


The equality in Eq. (E24) makes use of the fact that 


dpi 

dg 


E 


dA{{p[}) dp'i 

dp'i dg ' 


(E25) 


~Wf ~ P'iiddi — {Hc)pf) as derived in Eq. (D26). Evaluated at 
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g = 0, implies that p' = Pi, and therefore this gives 


dA 

dg 


= anBa, where 


g=0 

B„ = 


a 1-a 

2^P^d^ 


1-a 


((4)/3. - E,) . 


(E26) 

(E27) 


Recall that p^, qi are probabilties of the thermal states of H^, at inverse temperatures (3^, Ph respectively. With this, we can write 
the expansion of A with respect to g as 

A = l + angBa, + Q{g^). (E28) 


Later on, we will also need to evaluate the derivative of B^ w.r.t. a. This quantity, when evaluated at a = 1, has a close 
relation to the change in average energy of the cold bath (per copy), . 

Lemma 10. Let 


AC'(/3e) := 


3=0 


where recall j5f = f3c — g. Then 


B[ = 


dBo 


da 


Pc - Ph 


AC7'(/3/) = (/3e -/3/*) • var(4k. 


Ql^l 


Proof. From the definition of AC (Eq. (C14)) and using Eqs. (E18), (E19), (E21), we have 

AC 

-= tr[('r/3^ - TfsjHc]. 


Recalling that Pf = Pc — g and using Eq. (D36), from Eq. (E31) it follows 


(E29) 


(E30) 


(E31) 


-AC'(/3e) = -^ =--^ 

n n dg n dPf 

Now, let us evaluate the partial derivative of Ba w.r.t. a. Denoting 
Eq. (E27) 


/3/=/3c 


P± 

1i ’ 


= var(iJe)/3,. (E32) 

and invoking the chain rule of derivatives for 



. i 


^ q^rf Inr, ((Hc)/3, - E^ 

_ i 

. i 





(E33) 

(E34) 


Substituting a = 1 into Eq. (E33), we obtain that “ = 1- Also, J^iPfdl “((^c)/ 3 c — Ei) = 0 while the factor 

multiplied in front is finite. Therefore, we are left with the terms 


B'l = ^pjlnr* - E, 

i 


In ^ + iPh - Pc)Ei 


{{He) 


= {Pc - Ph)yar{He)p^ 
= ^E^^J\C'{Pc). 


Ei) 


(E35) 

(E36) 

(E37) 

(E38) 


The second equality comes from substituting = ^ = Pd)Ei . z^/Zc- In the third equality. In ^ is brought out 
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of the summation, while the summation yields 0. Subsequently, we invoke ~ Ei) = = 

-var(Hc)/3e- n 


ii) The expansion ofWa in the quasi-static heat engine 


We now proceed to derive an expansion of Wa valid for small g, and e. Note that Wi is defined through continuity to be the 

limit of the Renyi divergences at a —> 1, and the small £ and g expansion does not hold for a = 0, we shall have to examine Wi 

and Wq separately. 

(A) For e > 0, a G (0,1) U (1, oo) 

Wa = -s-j ^—— [ln(A-£“)-aln(l-e)] (E39) 

Ph\C^ Ij 

= -irj- —7T [in (1 + angBa + 0(g^) -£“)-« ln(l - e)] (E40) 

= [ctngBa + 0(g^) - £“ + 0(£^“) + Q{ge°‘) - a (-e + 0 (e^))] , (E41) 

= TpT-pY [oingBa — £“ + cue] + &{g^) + 0(e^“) + 0(pe“) + 0(e^). (E42) 

Ph[o: — Ij 

In the second equality, we have used the expansion of A derived in Eq. (E28). In the third equality, we use the Mercator series 


°° f-life+i 

ln(l + a;) = ^- - -|a;| < 1, 

k=l 


(E43) 


to expand both of the natural logarithms in line Eq. (E39). The order terms of Q{g^), Q{g^), Q{g^e°‘) vanish because they are 
of higher order compared with 0(p^) and Q{ge°‘). The last equality occurs because cQ{g{x)) = Q{g{x)) for any c G ]R\0. 

(B) For e > 0, a = 1 

We are now interested in finding a small e > 0, g > 0 expansion for Wi, which is defined through continuity of the Renyi 
divergences. Going back to Eq. (E12), note that Wi is the maximum value such that Eq. (E12) holds with equality, when all Da 
terms in Eq. (E12) are evaluated at a —?► 1. Recall that limQ,_j.i Da{p\\T) = D{p\\t), the relative entropy we have derived in 
Section D (see Eq. (BIO)). Therefore, one can write an equation for Wi in a more compact form: Wi is the value such that 






+ (1 - £)Wi - ^h 2 (e), 


(E44) 


where {Hc)p^ is the mean energy evaluated at temperature Tcoid, S{Pc) is the von Neumann entropy of the state and h 2 (£) 
is the binary entropy function. Rearranging Eq. (E44), we get 

Wi = ^ 

We can expand (E45) using a power law expansion in g and ^ for the terms in Eq.(E45), obtaining 


n(iJe)/3, - n(4)/3, - n— (^(/3,) - 5(/3/)) + —h2(£) 


Ph 


Ph 


(E45) 


Wi = [1 + £ + 0 (£ 2 )] 




n- 


dg 


5 + 0(5^) +Ah2(£) 

g =0 Ph 


(E46) 


To proceed, we recall that Pj = Pc — 9 and evaluate the term 


d{-{Hc)0,+pj:^sipf)) 


dg 


di{Hc)p, - PpSiPf)) 


Pc - Ph 


dpf 

var{Hc)0^- 


= -var(iTc)/3<, + ^vaT{Hc)p, 


Pf—Pc 


(E47) 


3=0 


Ph 


(E48) 
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This implies that when fully expanded, Eq. (E46) reads as 

=n-9 ^‘' a ^^ var(gc)/3, + /3;7^h2(e) + Q{eg) + 0(e)h2(e) + &ige'^) + 0(e^)h2(e) (E49) 

Ph 

+ 0 ( 5 ^) + (E50) 

=?T-5^^^^-^var(iTe)/3, + /?^^(-elne + e) + Q{eg) + 0(£^ Ine) + 0(e2) + 0(p^), (E51) 

Ph 

where we have used h 2 (e) = —elne + 0(e), which follows from finding the power-law expansion of the second term in Eq. 
(AlO). 


Although Eq. (E42) is not defined for a = 1, we can evaluate it in the limit a —1 to see if it coincides with the correct 
expression of VEi (in Eq. (E51)) at least for the leading order term (found in square brackets of Eq. (E42)). Eor the leading order 
term of Eq. (E42), we find 


lim — - - 

j3h[a - 


1 ) 


[angBa 


e + as] = 






aBa e — ae 

ng lim-lim - 

a—^l Ct — 1 a—>1 Q! — 1 

aBn 


ng lim -- + (—elne + e) 

a—>-l a — \ 

= ng ^‘' n ^'' var(gc)/3„ + ^^^(-elne + e). 
Ph 


(E52) 

(E53) 

(E54) 


The last equality holds because 


OiBq^ 

iim -- = lim —— 

oi^i a — 1 da 

= Wc - Ph) ■ var(Hc)/3,, 

where Eq. (E55) is derived from L’Hospital rule {Bi = 0 follows from the definition, see Eq. (E26)), and Eq. (E56) comes by 
invoking Lemma 10. Thus noting that Eq. (E54) is simply the first two terms in Eq. (E56), we conclude that the small g > 0 
and e > 0 expansion of Wa for a > 0 can be summarized as 


(E55) 

(E56) 


= 

f fa(a-i) l<^^9Sa - + ae] + Q{g^) -f 0(e2“) + 0(5fe“) + 0(e2) if a > 0, a 7^ 1 

I lim —— [anpi^o, — e“ + ae] + 0 (ep) + 0 (e^ Ine) + 0 (e^)-I- 0 (p^) ifa = l. 

I a^i+ Bhia - 1 


(E57) 

(E58) 


(C) Eor a = 0 

We will now investigate the a = 0 case. This is also particularly important to understand the difference between perfect 
and near perfect work, since in Section El, the impossibility of extracting perfect work arises from evaluating the allowed 
values of Wext under the a = 0 constraint. We show that by allowing e > 0, Wext > 0 is allowed once again. Recall 
E'o(p||<?) = Da{p\\q) = Thus fromEq. (E12) 

T>o(Pwlkw) - T'o(Pwlkw) > nDo(T/3^l|T/3j - nDo{TpP\Tp^) = 0. (E59) 

where the last equality follows from the fact that thermal states have full rank. This inequality is satisfied for any value of 
Ws^i, since whenever e > 0, is a full rank state, and £)o(PwIIav) = 0- Eurthermore, £)o(PwliAv) > 0 because all Renyi 
divergences are non-negative. Therefore, taking into account Eqs. (E57) and (E59), for quasi-static heat engines which extract 
near perfect work, we only need to solve 


kEext = inf W^, 

a>0 


(E60) 


where Wa is given by Eq. (E57). 
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d. The choice of e determines the infimum to evaluating Wtxt 


In this section, we will show that the infimum over all a > 0 in Eq. (E60) can be simplified to taking the infimum over a > R 
instead, where the parameter R determines how quickly e goes to 0 w.r.t. the parameter g. We define k in Lemma 11 and show 
its existence, for any function of e{g) such that limg_j,Q+ e(p) = 0. 


Lemma 11. For every continuous function e{g) > 0 satisfying limg_^Q+ e{g) = 0, 3 k S K>o s.t. 


^ I 0 ifK>K 

S{k) = lim - — = < (7 > 0 if K = R 


3-S-0+ g 


(E61) 


OO if K < K 


where n = +oo is allowed (that is to say, limg_j,o+ ^ diverges for every k) and a = +oo is also allowed. 


Proof The main idea in this proof is to divide the non-negative real line into an infinite sequence of intervals in an iterative 
process. We specify the ends of these intervals by constructing a sequence and evaluating 5 at these points. We then 

prove that according to our construction, there are only two possibilities; 

1) Ki forms a convergent sequence, where the limit lim„_).oo rin = R, or 

2) the ends of these intervals extend to infinity. In this case, R = oo. The way to construct this interval is as follows: in the 

first round, pick some ki > 0. The corresponding interval is [0, ki\. Evaluate If <5(/ci) = oo, then proceed to look at the 

interval [ki, Otherwise if di^Ki) < oo, choose K 2 = ^ and evaluate S{k 2 ). Depending on whether 6 {k 2 ) goes to infinity, 
we pick one of the intervals [0, K 2 ] or [k. 2 , ki]. 

A general expression of choosing can be written: during the nth round, define the sets , si°°^ such that 

5^°^ = {Ki\l < i < n and 6{Ki) = 0} 

< i <n and 6(^Ki) = 00 }. 

Note that if we find S{Ki) = c ^ 0 for some finite constant c, then our job is finished, i.e. R = Ki (We prove this later). 
Subsequently, define for n > 1, 


= min k 


k,GS. 


(0) 


and 


= max K. 


kGS- 


(oo) 


If either sets are empty, we use the convention that the corresponding minimization/maximization equals 0. Once these quantities 
are defined, we can choose the next interval by evaluating 




(E62) 


In the n-th round, the corresponding interval is Kn+i]. 

Let us now analyze why we can use this scheme to find R. Eirstly, consider the case where S{Ki) whenever evaluated, produces 
infinity. This means that in each round, increases with n (by the iterative scheme), and Kn^ = 0 always stays at zero. 

Note that this scheme has been constructed in a way such that lim„_>oo Kn = 00. Indeed, for all n, by using Eq. (E62), 


3 

^n+1 — 2^^ — 



(E63) 


which tends to infinity as n goes to infinity, whenever ki > 0. Later we will prove a property of the function 6 , which combined 
with this scenario means that 6 {k) = 00 for every k > 0. Therefore, R = 00 . 

Next, suppose that there exist an n-th round, such that = 00 and (5(n„+i) < 00 , as illustrated in Fig 6. Note that the 

function 6 {k) has a peculiar property, i.e. we know that if (5 (k„) = 00 , then for any k < Kn, 


S{k) 


lim (p) 


£”"(g) 

9 


= -foo. 


(E64) 


—>■+00 


—>•00 
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S(K) 


8„=oo 

-« 

for all K < K„, i 

6(k)=oo i 


s 


0 - 


T 8„+i<t» 


8n+2 ? 


I for all K > K„+|, 

i 8( k )=0 

-i-1-»-K 


FIG. 6. Illustration of the scenario where S(Kn) = oo and S(k„+i) < oo. 


On the other hand, if S(Kn+i) 


0, then we know that for any k > /c„+i, 

S(k) = lim £«-""+! (5) =-^ = 0. 

-J-O 


(E65) 


Moreover, if S(Kj) = c ^ 0 for some positive, finite c, then following the same arguments, one can easily see that for all k < kj, 
S{k) = 00 and for k > Kj, 6 {k) = 0. In this case we find that R = kj. These observations are illustrated in Figure 6 for clarity. 

One can now evaluate «;„+2 (which is the midpoint of and and its corresponding value of (5 (k„+ 2). From this point 

on, in each iteration we either find R exactly (whenever the function 6 when evaluated produces a finite, non-zero number), 
or the length of the next interval gets halved, and goes to zero in the limit of n —> c». This, by Eq. (E62), also implies that 
lim„_>oo = lim„_>,oo We also know the following: 

1) for all K < k^\S{k) = 00 , 

2) for all K > Kn \S{K) = 0. 


Therefore, we see that k exists and k = lim„ 




= lim„ 


JO) 


. By this we conclude the proof. 


□ 


To provide some intuition about how R compares the rate of convergence e, g —> 0, let us look at the following examples: 

1) Consider ei(p) = exp(—1/p). Then R = 0 with a = 00. 

2) Consider £2(5) = plnp. Then R = \ with a = 00 . 

3) Consider £3(5) = c • for fc > 0. Then R = k with cr = c. 

In the next lemma, we consider the scenario of near perfect work, given in Def. 2, and show that this imposes a finite range of 
values R should take. Given a particular R, we also show that the minimization of Eq. (E60) changes with R. 

Lemma 12. Given any e{g) £ (0,1] as a continuous function of g, where p > 0. //limg_^Q+ £(p) = 0 and limg_^Q+ = 0, 
then the following holds: 

1 . The quantity R (defined in Lemma 11 ) can only have any value in R £ [0,1], where limg_j,Q+ = q iq hold if 

R=l. ® 

2. The extractable work can be written as 


, _ naBa 

inf -- 

ol>k q: — i 


fig) 


Wext = P • 

where limg_j,o+ /(p) = 0 and infQ,>K can be exchanged for infQ,>K ifR = 0. 

Proof Eirstly, let us use Eq. (E57) to simplify our expression for Wext: WJxt = infa>o ^a, where 

gWa + 0(p^) + 0(£^“) + 0(p£“) + 0(£^) if a £ (0,1) U (1, 00) 


(E66) 


PhWc. = 


gWi + 0(£p) + 0(£^ ln£) -F 0(£^) -f 0(p^) if a = 1, 


(E67) 
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and 


and for a = 1 


If e 

Wa :=- anBa + a- 

a - 1 V g 



Wi = 



anBa \ 
a-1) 



(E68) 


(E69) 


Erom now on, the order terms in Eq. (E67) can be neglected, since it can be checked that all of them are of higher order compared 
to the terms we grouped in Wa, in the limit of vanishing g. Even then, we note that due to the complicated form of VEext, it is 
not straightforward to begin our proof with the assumption limg_^Q+ = 0. 

Instead, we begin by noting that given a function e{g) that satishes the conditions of the above lemma, then one can invoke 
Lemma 11, and therefore there exists aRG ffi>o such that Eq. (E61) holds. We then, for all possible k G K>o, evaluate all TEq- 
to take the infimum and obtain lEext- Given lEext, we then evaluate the quantity limg_j,o+ = 0. 

The value of R determines how the limits of quantities like |^ behave. Therefore, we need to split the analysis into three 
different regimes: R G [0, 1 ), R — 1 , R G (1, oo). 

1) Eor R G [0,1) 

Eor this case, we know the following limits: 

A. limg^o+ 1 = 0- 

B. Eor a < R, limg_j,Q+ ^ = oo. 

C. Eor a = R, limg_j,Q+ ^ = a > 0 . 

D. Eor a > R, limg_j,o+ ^ = 0. 

E. Note that 3 ki > R such that 1 — fci > 0. Thus limg_^Q+ | Ine = limg_j.o+ ^ Ine = 0 

if a G [0, R) 
ii a = R 

if a G {R, 1) (E70) 

if a G ( 1 , oo) 
if a = 1, 


Therefore, by using Eq. (E68) and (E69) (for a = 1 separately) we have 


Wa = 


+00 

anBa 
a — 1 
anBa 
a — 1 
anBa 


a —I 
anBa 

lim 


0 


0(^ 

9 

0(^ 


o!->i a — 1 


9 

+ 0 


elne 


where the expression in Eq.(E70) has been written as a leading order term, plus higher order terms that vanish in the limit 
9^0. 

Therefore, we conclude that for R G [0,1) and any cr > 0, due to continuity in a of , 


/3/iWext = ini Wa = 9- 

a>0 


. _ anBa 

ini - 

cx>K a — 1 


+ 0 (/( 5 )) 


(E71) 


where / satisfies limg_>o+ /(p) = 0 in the expression of Eq. (E70), Both functions vanish as g tends to zero. Note that if At = 0, 
then inf„>g can be exchanged for infa>R since in Eq. (E60) the point a = 0 was already excluded. 

We can now calculate limg_^o+ for R G [0,1) and any ct > 0: 


AS — elne — (1 — e) ln(l — e) 

lim - = lim - -z -r- 

9^0+ fEext 9^0+ g 


lim - 

9-A0+ infa>s 


OLTlBfy 

a-1 


( 

eine 

\^0 (Item E) 


\ 

e + 0(£^) 

9 

—f 0 (Item A) / 


(E72) 
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where we have assumed that 

i„t 25% > 0, 

a>K, Q — 1 

As we will see later (see Eq. (E89)), Eq. (E73) holds if k > 0. However, 

anBa 


a — 1 


if a = 0 and we need to use Eq. (E70) for the case a € 1) for k = 0. Erom which we conclude that 

/3/tWext = inf Wa > = 0(/(g)), 

a>0 


(E73) 


(E74) 


(E75) 


thus we have 


AS' 

lim - < 

g—>0+ lEext 


lim 

g-*' 0 + 


—elne — (1 — e) ln(l — e) 

mTJJ 


lim g(—_ £i/2t _ 

S-*'0+ 


(E76) 


thus from Eqs. (E74), (E75), and (E76), we conclude that Eq. Eqs. (E7 1 ) and (E72) are still valid when k = 0. To summarize, so 
far we have proven that whenever k G [0,1), Eq. (E66) holds for some /(g) which vanishes as g tends to zero, and furthermore 

C = 0 . 

2) Eor K € (1, oo) 

In this regime, like the previous analysis, we can list out the following limits: 

A. limg_,o+ 1 = 0. 

B. By definition of R, for a < 1, limg_^o+ ^ = oo- 

C. limg_^Q+ = oo since both | and Ine goes to infinity as g —>■ 0. 

Therefore, by using Eq. (E68) and (E69) (for a = \ separately) we have 

f i + ®(5)] if « e 

= I f • [-£ln£ + 0(e) + 0(g)] if a = 

[ae + 0 (e“)+ 0 ( 5 )] if « e 


[ 0 , 1 ) 

1 (E77) 

(1,00). 


Note that for all of these expressions in Eq. (E77), Wa 00 . Next we want to calculate Wext, which is the infimum of Wa, 
taken over all a > 0. Note that in the limit of vanishing g, e also goes to zero. Therefore in Eq. (E77), the equation of Wa which 
vanishes most quickly in the limit g —> 0 happens when a G (1, 00 ). Therefore, we conclude that for R G (1, 00 ) and any cr > 0, 


/^hWext = inf lEa = g • 

a>l 


inf - 

a>l a — 


£ 

15 


+ 0 (/(<?)) 


= e + g • 0(/(g)) 


(E78) 


We can now calculate limg_).g+ ^ for k G (1, 00 ) and any cr > 0: 


lim 

9-r0+ 


W 


lim 

g-rO+ 


—elne — (1 — e) ln(l — e) 
£ 


£ln£ £ + 0(£^) 
hm- ^—- = +03. 

y _£ 

—>•00 —>-1 


(E79) 


Erom this, we note that the whole regime of k G (1, 00 ) does not contain any cases corresponding to our condition of interest: 
limg_).o+ = 0 never holds. 

3) Eor R = 1 

Similar to the first two cases, we again list out the relevant limits: 

A. limg_j,Q+ 1 = 0 ’ for some cr > 0. 

B. Eor a < 1, limg_j.o+ ^ = 00 . 

C. Eor a > 1, limg_j.o+ ^ = 0. 














Therefore, by using Eq. (E68) and (E69) (for a = 1 separately) we have 
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i • [-eln£: + 0(e) + 0(p)] 


if a G [0,1) 

if a = 1 && <7 > 0 


W^ = 


lim 


aBn 


eln£ aBr, 

- > n lim 


1 

a-l 


if a = 1 && (7 = 0 
if a G (1, oo). 


(E80) 


*1 a — 1 g a-fi a — 1 

anBa + aa ~ 0 (V) 

Note that for a G [0,1) and the case a = 1 && a > 0, Wa tends to infinity, while for the other cases Wa is finite. 
Therefore, we can conclude that for S = 1, 


/3hWext = g ■ 


where f{g) = — vanishes as g tends to zero. 


inf -- 

■Dt>l Q( — 1 


{nBa + cr) + 0 ifig)) 


Now, we evaluate the limit limg_^Q+ ^ for k = 1 and any tr > 0; 

AS —sins — (1 — e) ln(l — e) —sins s + 0(s^) 

lim - = lim . ---r—- = lim- 

g^o+ W 9^0+ |^inf„>i ^ (anF^ + aa)j g 3^°^ ^ ' 9 , ^'9 , 


(E81) 


(E82) 


This limit of interest can be zero if and only if limg_j.Q+ = 0. 


We have calculated the limits limg_j.o+ AS/W^a to leading order in g for all functions e{g) > 0 satisfying limg_j.g+ s = 0. 
These are found in Eqs. (E72), (E79), and (E82). We have found that limg_j,g+ AS/W^xt = 0 occurs only in two cases: 

i) K G [0,1), and 

ii) K = 1 and limg_j.o+ = 0. 

The amount of work, Wext is found in Eq. (E71) and (E81) respectively. Indeed, they take the form of Eq. (E66), for different 
functions f{g). With this, we conclude the proof of the lemma. □ 


e. Solving the infimum for Wext 


We have seen in Lemma 12 that the function corresponds to the largest order term in Wext w.r.t. small g (quasi-static heat 
engine). Our next objective is to find the infimum of over a G [k, c»] appearing in Eq. (E66). Such an infimum is is not 
easy to evaluate, but whenever the cold bath consists of multiple identical qubits, we show that the derivative has some 

nice properties. Roughly speaking, we show that this derivative does not have many roots, which in turn means that does 
not have many turning points. We have used this to prove in Lemma 14 that the infimum is either obtained at a = it or a —> oo. 


The derivative of w.r.t. a is given by 

d aB„ Bn B'„ 
-I- a 


aBn 


BL 


da a — 1 a — 1 


— 1 {a — 1)^ {a — 1)^ 


Bn 




BL 


{a - 1): 


;G{a), 


where 


Bn 


G{a) :=a{a-l)- —. 


(E83) 


(E84) 


Now, we shall evaluate the quantities Ba,Bn, and for the case of qubits (see Assumption (A.5)), where the energy 
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levels are {0, E}. By using Eq. (E18), we evaluate the quantity Ba defined by Eq. (E27) to obtain a simple expression: 


g-/3c£; 

1 pOiPhE 

_ ^ _ _ _ _ 

1 -|_ g/3c-E ^aPhE _|_ g(/3h+a/3c)-E 

E \ e“^'*'®(l + e'^‘=-®) 

X _|_ g/3c-E ^ ^aPhE _|_ g(/3fe+a/3c)-E 
g(/3h+a/3c)-E _ g(/3e+a^fe)£; 

X _|_ g/3c-E ga/3h,-E _j_ g(/3h+a/3c)-B 


(E85) 

(E 86 ) 

(E87) 

(E 88 ) 


We note that Eq. (E 88 ) is zero only if a = 1, and thus for a ^ 1, aBa/{a — 1) 7 ^ 0. Erom Eq. (E56), we know that 
limQ_>i aBal{a — 1 ) > 0 , thus due to continuity, 


aBa 
a — 1 


> 0 


V a > 0. 


(E89) 


We also derive the first derivative of B^ w.r.t. a for the special case of qubits: 

dBa, E^iP, - ph) 


BL = 


{Ph+aPo+aPh)E 


dOL \^^aPhE _j_ g(/ 3 h+a/ 3 c)iJj ^ 


(E90) 


Note that since /3c > Ph by definition, therefore whenever E > 0, then B'^ > 0 always holds. By further algebraic manipulation, 
we compute the first derivative of the function 


_ cosh[w{Pc,Ph,a)E] 
da B'^ cosh{PcE/2) 


(E91) 


where w{Pc, Ph, a) = {Pc - Ph)a + Ph- 

We have written Eq. (E83) in this form, since for the special case of qubits, namely Eq. (E90), B'^ > 0 is always true. 
Therefore, looking at the function G{a) whether it is positive or negative) will tell us whether (and therefore Wa) is 
increasing or decreasing in a particular interval. 

In Lemma 13, we identify the conditions on the energy spacing E such that several different properties of G{a) hold. 


Lemma 13. Consider G{a) = a{a — 1) — where Ba, B{^ is defined in Eq. (E 88 ) and (E90). Then the following holds: 
1) IfE{P, - Ph) tanh(/3cE;/2) > 2, 


3 0 < T < 1 s.t. 


G{a) <0 Va S (r, 1) U (1, 00 ) 


2) IfE{p, - Ph) ta.nh{P,E/2) < 2, 


3 a > 1 s.t. 


G(a)>0 VaG (0, l)U(l,a) 
G{a) <0 Va G {a, 00 ). 


(E92) 


(E93) 


3) IfEiP, - Ph) ta.nh{P,E/2) = 2, 


G{a) >0 Va G (0,1) 

G(a)<0 VaG(l,cx)). (E94) 


Proof. Eirst we note that since Bi = 0, therefore G(l) = 0. Let us also compute the derivative of G(a) w.r.t. a: 

cosh {{-pc/2 + Ph + iPc - Ph)a)E) 


G'ia) = 2a - 1 - 


cosh{PcE/2) 


(E95) 


Before we continue, there are several properties of the function G'(a) which we shall make use of. Eirstly, note that G'(l) = 0, 
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in other words, G' has a root at a = 1. Also, G'{oo) = —oo for any value of E' > 0, Ph > 0, /3c > Ph *■ Also, since 2a — 1 
is linear (and hence both convex and concave), while the — cosh function is strictly concave ^ , therefore the function G'(a) is 

strictly concave. This implies that the second derivative G"{a) = is strictly decreasing w.r.t. a. 

The properties of G'{a) indicate that we can fully analyze the function by considering 3 different cases: 

1. G' has two roots at a = {a, 1}, wherewhere a G (—oo, !)■ This corresponds to the case G"(l) < 0. 

2. G' has two roots at a = {1, a}, where a G (1, oo). This corresponds to the case G"{1) > 0. 

3. G' has a single root at a = 1. This corresponds to the case G"{1) = 0. 





FIG. 7. G"(l) < 0. 


FIG. 8. G"(l) > 0. 


FIG. 9. G"(l) = 0. 


FIG. 10. A convex function G'{a) and its corresponding G{a), for different values of G"(a). 


We shall now consider these cases one by one. Suppose that 


G"(l) 


G"(a) 




2 - (/3c - /3h)Etanh < 0, 


(E96) 


then G''(a) < 0 for all a G (1, oo). Note that Eq. (E96) corresponds to the first condition in the lemma stated above. 

This information about the second derivative G''{a) now allows us to conclude the following about G(a): 

1. If for all a G (1, oo), G''(a) < 0, then we know that G'(a) < 0 holds for all a G (1, oo) too. Eurthermore, this implies 
that G(a) is monotonically decreasing in the interval (1, oo) and therefore, G(a) < 0 for all a G (1, oo). 

2. G"(l) < 0 also implies that there exists an interval (r, 1) such that G'(l) > 0 (See Eig. 7). And since G(l) = 0, this 
implies that within the interval (r, 1), G(a) < 0. 


With this, we prove the first statement of the lemma. 

Let us now analyze the second case, where G"(l) > 0. This implies that G'(a) > 0 at least for some interval a G (1, a), then 
G'(a) changes sign exactly once at a = a, and goes to —oo. (Refer to Eig. 8). Also, recall that in the limit of a —oo, G also 
goes to —oo. Therefore, we conclude that there exists some a such that 


G(a) 


> 0 a G (1, a) 
< 0 a G (a, oo) 


(E97) 


With this, we prove the second statement of the lemma. 

Einally, we look at the case where G"(l) = 0, and make the following observations: 

1. Since the function G'(q;) is concave, and since G"(l) = 0 implies that a = 1 is an extremum point for the function G'(a), 
we know that it must also be the global maximum. Therefore, we know that for any a ^ l,G'(a) < 0. 


* This is due to the fact that 2a increases linearly w.r.t. a, while the cosh term increases exponentially. 

^ To be more precise; due to the concavity of f{x) = —acosh(f) + xc) for a > 0. This follows from the strict concavity of the cosh function, the invaiiancy 


of strict cone 
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2. Since for the interval a S (— cxd, 1), G'{a) < 0 and we know that G(l) = 0, therefore we can deduce that for any 
a S (—c», 1), G{a) > 0. 

3. Since for the interval a G (1, oo), G'(a) < 0 and we know that G(l) = 0, therefore we can deduce that for any 
a G (l,oo), G(a) < 0. 

With this, we prove the final statement of the lemma, and complete our proof □ 

To summarize, in Lemma 13 we have identified conditions involving the energy gap of H^, and the temperatures fin, Pc- 
Depending on whether these conditions are satisfied, we can describe the positivity/negativity of G(a) for different regimes of 
a. Comparing these different scenarios, we prove in Lemma 14 that for a quasi-static heat engine, the minimum of infQ,>K 
is obtained only either at a = k or a = oo. 

Lemma 14. There exists some Q < v < 1 such that V k satisfying v < n < 1, the following infimum is obtained at one of two 
points 


inf 


aBn 


y->K a — \ 


= inf < lim 


aBc 


, lim 


olBc 


*K a — 1 a—^OO Q; — 1 


^ aBa 

< lim - 

a—>/3 a — 1 


y p G { k , oo), 


where B^ is defined in Eq. (E88). Furthermore, if E{Pc — Ph) tanh(/3ci?/2) < 2, then we can set z/ = 0. 

Proof. 1. If . 

d aBa j > 0 V q; G (0,1) U (1, a) for some a > 1 

daa — 1 1 < 0 V a G (a, oo). 


then V K G (0,1), 
inf 


aBa 


a>fi a — 1 


= inf < lim 


aBa 


lim 


aBa 


a^K a — 1 ' Q—>oo a — 1 I a—>/3 a — 1 


< lim 


aBa 


V /3 G (k, oo). 


(E98) 


(E99) 


(ElOO) 


Recall from Eq. (E83) that 


BL 


d aBa 


da a — 1 {a — 1) 


rG(a), 


(ElOl) 


where B'a > 0, and we have derived some properties of G{a) in Lemma 13. In this proof, we apply Lemma 13 directly to 
consider the three scenarios detailed in Lemma 13. 

Eirst, consider the first statement of Lemma 13. If E{Pc — Ph) tanh(/3cG/2) > 2, then 3 0 < f < 1 s.t. 


d —a Q y Q, g (f, 1) U (l,oo) 

aa a — 1 


(E102) 


then by continuity of in a, we conclude that V k satisfying t < k < 1 


. r fxBa ,. aBa 

ini -= lim - 

a>K a — 1 ct—>oo a — \ 


. p I ,. aBa ,. aBa 

= mt < lim -, lim - 

->■« a — \ a —>-00 a — 1 


< lim y P G {k,oo). 

a^/3 a — I 


Next, consider the second and third statements of Lemma 13 jointly, where E{Pa — Ph) tanh(/3c£’/2) < 2. 
statements proved in Lemma 13 (namely, Eq. (E93) and (E94)) can be rewritten as the fact that there exists a > 


(E103) 

Note that both 
1 s.t. 


d aBa 1>0 for a G (0,1) U (1, a) 
da a — 1 I < 0 for a G (a, oo). 


(E104) 


In fact, the third statement is simply a special case of the second, where a = 1. If Eq. (E104) holds, then V k G (0,1), 

inf = inf ( lim lim < Hm V /3 G (k, oo). (E105) 

a>K a — 1 l^ct—a — 1 Q—>oo a — 1 j a— )-/3 a — 1 


By setting r = 0, we see that the statement of Lemma 14 is achieved. 

Therefore, since we have analyzed all three cases stated in Lemma 13, we conclude that there always exists v G [0,1) such 
that Eq. (E98) will always be satisfied V k G {v, 1). 

□ 
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f. Main results: evaluating the efficiency 

In this section, we derive the efficiency of quasi-static heat engines in the nano /quantum regime. We first need to define the 
quantity 


n ■ = 


min 


EdPc-Ph) 

} 1 -f ’ 


(E106) 


where recall that Ei is the energy gap of the cold bath qubits, as described in Eq (E18) and the sentence right after it. Recall that 
n denotes the number of qubits in the cold bath, where n G is any positive integer. Before stating the maximum efficiency, 
we will derive the efficiency as a function of R, defined in Lemma 11 (recall that this parameter is determined by the choice of 
e). Eor simplicity, we will still consider the special case where Ei = E for all i in Lemma 15, (i.e. all qubits of the cold bath are 
identical). Lemma 15 shows us that under the condition of extracting near perfect work, one can choose e (and therefore R) such 
that a certain maximum efficiency value is achieved. The closer R is to unity, the slower limg_>o+ /S.S/W converges to zero, and 
also the closer the efficiency is to the Carnot efficiency. 

Using this lemma, we prove the achievability of the Carnot efficiency which depends on 17. This is the main result of our 
work, which is stated in Theorem 2. 

Lemma 15 (Quasi-static efficiencies as a function of R). For any n G Z’*' number of qubits, consider quasi-static heat engines 
(Def 3) as a function ofR (defined in Lemma 11) which extract near perfect work (Def. 2). For any k G (0, oo)\{1}, define 

j(k) := ^ (E107) 

K — 1 

where is defined in Eq. (E 88 ), while 7 ( 1 ) and 7 (cx)) are defined by taking the limits k —> 1, cx) respectively. 

If Ft < 1 (see Eq. (E106)).- 

1) There exists v G [0,1) such that for any R G (v, 1] (and\im.g^Q+(e\B.e) / g = 0 ifR = 1), the maximum efficiency is 

= 1 + o + 0(/(ff)) + 0(5) + 0(£), (E108) 

Pc - Ph 7(«) 

where 7 ( 1 ) > 7 (k) with equality iffR = 1 and limg_>o+ f{g) = 0 . 

2) The corresponding amount of work extracted is 

TL 

W,,fR)=g—Yi(R) + Q{f(g))]. (E109) 


7/17 > 1.- 

1) There exists v' G [0,1) such that for any R G {E, 1] (and limg_>o+ {elns)/g = 0ifR = 1), the maximum efficiency is 


7(1) 


g (k) = 1 -f 


Pc - Ph 7(00) 


0(/(5)) +0(5) +0(e), 


(El 10) 


where 7 ( 1 ) < 7 (cx)). 

2) The corresponding amount of work extracted is 


lEext(K) = g^ [ 7 ( 00 ) -f 0 ifig))] 

Proof. Lirstly, let us begin by deriving the explicit form for 7 ( 1 ) and 7 ( 00 ): 

7(1) = lim —= lim -f = f, e\2 ^ 

a—>-1 Q! — 1 a—>1 (1 -|- 


where we have made use of the L’Hospital rule. Lor a —> 00 , since 

R pPhE _ „/3,;E -a{0a-i3h)E 

lim B„ = lim 


E 


a—foo ^ a— foo 1 -|- p^l^cE ^(ihE _j_ g Q,(/3c $h)E ^ _j_ f^^cE ’ 


(El 11) 


(El 12) 
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therefore we have 


7 ( 00 ) = lim ("1 + 3^= ^ (E113) 

\ a — 1J 1 + 

By Lemma 14, we know that the infimum of 7 (a) for a G \R, 00 ) and R G {v, 1] is either at a = k or a —?> 00 . Therefore, if we 
take the ratio of Eqs. (El 12) and (El 13) to be 


7(1) ^ - Ph) 

7 ( 00 ) 1 + 


(El 14) 


then 7 ( 00 ) > 7 ( 1 ) > 7 (k), therefore the inhmum of 7 (a) for a G [it, 00 ) and R G (i^, 1] has to be obtained at a = R. Taking 
this into account and using the condition which is equivalent to that of near perfect work in Eq. (All), we can use Lemma 14, 
to calculate the amount of work extracted; 


W4xt 


inf Wa = g- 

a>0 


inf ® 

a>K 


= [7(tj) + 0 (/(s))], 


(El 15) 


where limg_>o+ f{g) = 0. On the other hand, we can calculate AC, which is the change of average energy in the cold bath 
system, (recall this is done by Taylor expansion around g = 0) 

AC = n - {E)l) g + 0 (g^) = £^5 + © (ff") ■ (El 16) 


Using Eq. (C15), we have AW = (1 — e) Wext- The (inverse) efficiency, according to the definition (Cl), is thus 

. AC _ , , ^7(1)/ {13c - Ph)g + 0 (g^) 

^ ^ fEext ^ n-f{R)g/l3h + Qigf{g)) 

“ ‘ (A^ ^ + e(/(s)) + e(s) + e(e), 


(El 17) 
(El 18) 


where we have used limg_,,o+ f{g) = 0 which is proven in Lemma 12. We will now investigate the efficiency when 0 > 1 is 
satisfied. Using U > 1 and Eq. (El 14), we have that 7 ( 00 ) < 7 ( 1 ). Thus from Lemma 14, due to continuity in R of 'y{R) we 
conclude that there exists av' G [0,1) such that for any R G {y', 1], 


inf 7 (a) = 7 ( 00 ). (E119) 

CX.>K 


Therefore, since we are considering near perfect work, Eq. (All) holds and we can use Lemma 12 to calculate the amount of 
work extracted 


Wext = inf Wa = p ■ 

Q ;>0 


inf 

“>« Ph 


= 9 


Ph 


7 ( 00 ) + —fig) 


where limg_>o+ fig) = 0. Thus using the definition of inverse efficiency (Eq. (Cl)), together with Eq. (El 16), we have 


^ AC _ njil)/il3c-Ph)g + &{g‘^) 

(Eext ^ njioo)g/l3h + Qigfig)) 

^ ^ (R R 1 \ ^^•I'^d)) + 0(ff) + ©(^)) 

iPc - Ph) 7 ( 00 ) 


where we have used limg_,,Q+ fig) =0 which is proven in Lemma 12. 


(E120) 


(E121) 

(E122) 

□ 


We will now use Lemma 15 to conclude our main result of this letter. 


Lemma 16. Consider the case of near perfect work (Def. (2)) and all cold bath qubits are identical (i.e. Ei = E for i = 
1 ,..., n), then: 
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1) < 1 (see Eq. (E106)j the optimal achievable efficiency r^max (see Eq. (C4)j is the Carnot efficiency: 

’'"■> = (‘ + 4 ^)" 

What is more, this efficiency is only achieved for quasi-static heat engines, i.e. rfmax = (see Eq. (C8)j. 

2) If Cl > \ and the heat engine is quasi-static, then the optimal achievable efficiency is (see Eq. (C8)j 

= '■ 

3) If Cl > 1 the maximum achievable efficiency rymax (see Eq. (C4)j, is strictly less that the Carnot efficiency. 

Proof. In Lemma 5, we found that the Carnot Efficiency is an upper bound for the efficiency when we are extracting near perfect 
work. We also found that Eq. (All) is satisfied iff we are extracting near perfect work. In Lemma 15, we derived the optimal 
achievable efficiency for quasi-static heat engines as a function of R when Eq. (All) is satished. By choosing k < 1 arbitrarily 
close to one, if (1 < 1 is satished, we will thus achieve an efficiency arbitrarily close to the Carnot efficiency. Thus since the 
upper bound is equal to the lower bound, we prove part 1) of the Theorem. Part 2) of the Theorem follows from setting k = 1 in 
Lemma 5 when 12 > 1 is satished. □ 

By making use of Lemma 15, one can generalize Lemma 16 to consider the more general case stated in A5 (at the begging of 
section E 2 c) where the cold bath still consists of qubits, however the energy gaps of the qubits can be arbitrary. Lor convenience, 
we re-write the general cold bath Hamiltonian here: for a set of variables Ei > 0, • • • , En > 0, 

n 

= where H^ = Ek\E){E\, (E125) 

fe=i 

Under the more general form of the cold bath Eq. (E125), we have the following theorem. 

Theorem 2. [Quantum/Nano heat engine efficiency] Consider the case of near perfect work (Def. (2)), when the cold bath 
consists of multiple qubits with energy gaps {Ei]^^^. 

1) If Cl < 1 (see Eq. (E106)) the optimal achievable efficiency rj^^x (see Eq. (C4)) is the Carnot efficiency: 

What is more, this efficiency is only achieved for quasi-static heat engines, i.e. p^ax = t?max (s^^ Eq. (C8)). 

2) If Cl > 1 and the heat engine is quasi-static, then the optimal achievable efficiency is (see Eq. (C8)) 

3} If Cl > 1 the maximum achievable efficiency rymax (see Eq. (C4)), is strictly less that the Carnot efficiency. 

4) Allowing for correlations between the final state of the battery and cold bath cannot improve the efficiencies achieved in 
1), 2) and 3) above. 

Proof. 1) is relatively simple to prove: as long as there exists a qubit with energy Ei such that < 1, one way to 

achieve Carnot efficiency is to simply disregard the rest of the cold bath, and act only on such qubits. The result is a simple 
application of 1) in Lemma 16. This strategy might not be optimal in terms of work extracted, but it is sufficient for our proof. 
Lor 2) and 3), suppose that U > 1. Since U is a monotonic function of E, we conclude that for all Ei where 1 < J < n, 
•= > 1- By Lemma 15, we see that this implies that the work extractable for all the individual qubits (which is 

an optimization problem over all a > 0) is obtained at a —oo. In general, considering the qubits collectively does not mean 
that the collective Wext is additive. This is because the minima of two functions is not necessarily the minima of these individual 
functions added together, as illustrated in the l.h.s. and middle diagrams of Ligure. 11. However, (as illustrated on r.h.s. diagram 
of Ligure. 11), when all the functions have their minima at the same value, then the collective minima is also obtained at that 
value. 
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FIG. 11. Illustration of the minima of two individual functions/(a:), g(a;) and minima of/(a;) + g{x). 


Next, we show that no matter which subset of qubits S one picks, Carnot efficiency cannot be achieved. We begin by 
introducing the notation 7 i(a), where 7 i{a) is defined similarly with 7 (a) in Eq. (E107) and (E 88 ), and the index i indicates 
that E is substituted by Ei in Eq. (E 88 ). Furthermore, recall that from Eq. (El 14), 17^ > 1 is equivalent to 7 ^( 1 ) > 7 i(c»). Now, 
consider any subset of qubit indices S, the amount of extractable work (as a function of g) is 




Ph 


^7i(oo)+/(p) , 

-iGS 


(E128) 


where limg^o+ fig) = 0 . 

On the other hand, we have that AC depends on the individual reduced qubit states, since there are no interaction terms in 
Hco\a- Therefore, similar to Eq. (El 16), 


AC-^ 


g 

Pc - Ph 


^7,(1) + 0(p2) . 

i^S 


(E129) 


Following the same proof in Eq. (E121) Lemma 15, 


?7 ^(k) = 1 + 


AC 


- s = 1 + 


Ph 

Pc- Ph E*G5 7i(oo) 


+ 0(p) + 0(£). 


(El 30) 


As we have observed before, the inverse of the Carnot efficiency 77 ^^ = 1 + - Furthermore, notice that by Eq. (El 14), the 

condition > 1 implies that 7 i(l) > 7 i(oo). Since > 1 is true for all 1 < * < n, therefore > 1- 

Lastly, part 4) is proven in Section F. □ 


Suppose n is large. Then since we have a spectrum which looks like a quasi-continuum: the full range of the spectrum is very 
large, compared to the individual energy gaps. One expects that in such a case, baths are of high temperature (small values of 
/3), then the effects of quantization should give us the classical observations of being able to achieve Carnot always. This can be 
seen, that for E^ia = min Ei, if the quantities PcEmim PhEmm ^ 1 , then 

iG{l, ' .n} 


T^min(/3c - Ph) 

^ —j— g /^c-fi'min 


^ EminiPc 


Ph) < 1. 


Whenever O < 1, we know that Carnot efficiency is achievable. 


(E131) 


3. Running the heat engine for many cycles quasi-statically 

We have so far proven that a heat engine can achieve the Carnot efficiency when O < 1. However, as like with macroscopic 
heat engines, this can only be achieved when the heat engine runs quasi-statically. Macroscopic heat engines can then extract a 
finite amount of work by running the heat engine over many cycles (in fact, over any infinite number of cycles if they want to 
obtain the Carnot efficiency in order to rvin quasi-statically). The following lemma, shows that when fl < 1, a nano-scale heat 
engine with a machine that runs over infinitely many cycles can also achieve the Carnot efficiency, while extracting any finite 
amount of work W with vanishing entropy increase in the battery. 
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A Top view 
of Fig. 8 


Hot bath 



Cold bath 

FIG. 12. Depiction (top view) of a heat engine comprising of a hot bath, 
a cold bath consisting of n identical qubits, a machine and a battery. In 
each cycle, the machine interacts specifically with one qubit from the 
cold bath, together with the hot bat and battery. After the end of one 
cycle, the machine is returned to its original state, and acts on a different 
qubit in the cold bath. 



FIG. 13. Side view of the heat engine. After each cycle of 
the machine, the battery, depicted here as a weight moves 
upward by a small amount. After N machine cycles, it 

has been lifted from its original position j Ej ^ to a final 
state that has most of its weight on \Ek)- 


For simplicity, we will work with the case in which the quasi-continuum battery has a part of its spectrum equal to that of at 
least N qubits, each with an energy gap Wext- We work within this subspace. We will run a heat engine between a hot bath, cold 
bath using a machine which performs N cyclic cycles. Let Ej and be the smallest and largest energy eigenvalues within this 
subspace respectively. We let the initial state of the battery be 

p'i = \E,){E,\, (E132) 

H^\Ej) = Ej \Ej) while we wish the final state of the battery to be of the form 

= r\Ek){Ek\ + {I - r) p^, (E133) 

where H^ff\Ei;) = p^ is some orthogonal state to \Ek) and the value of the probability r is to be specified in the 

following lemma. We will define the amount of work extracted from the machine for N cycles 

W,y,-.= Ek-Ej. (E134) 

Eor simplicity, we will consider the case that the cold bath consists of n identical qubits with 17 < 1, and during each cycle the 
machine interacts with one qubit from the cold bath. The running of the heat engine is depicted in Eig. 12 and 13. 

Corollary 2. [Many quasi-static heat engine cycles] Let W be the finite amount of work we wish to extract. Then for allW > 0 
and (5 > 0 there exists an n identical qubit cold bath (with 17 < Ij and an N € number of machine cycles with n> N such 
that: 

1) rjc (L V ^ Vc — where the efficiency p is the efficiency per cycle and is defined by Eq. (Cl), and rjc = ^ — Ph/Pc is 
the Carnot efficiency, 

2 } ILeyc > IL - ,5, 


3) S{p%) = 0, S{p\f) < 5, and 















46 


4) r > I — 6. 

whats more, 6 ^ 0 as N ^ +oo. 

Proof. Since in the qubit subspace, the spectrum is that of at least N qubits, we can write the initial state in the form 

pl = \E,){Ef’^^, (E135) 

with iJw \Ef)'^^ = Ej \Ej)'^^. We can now apply the heat engine results of Lemma 15 to the setup. Namely, we can apply 
the results of one cycle to each of the qubit subspaces of the battery in parallel. From Lemma 15 we conclude that this can 
be achieved with an efficiency given by Eq. (E108) and extract an amount of work per qubit/cycle given by Eq. (E109). For 
simplicity, we will run the heat engine using one qubit of the cold bath at a time. The hnal state of the battery is thus 

p\, = [(1 - e) \Eu){Ek\ + e \E,){Eff^ . (E136) 

Noting that \Ek){Ek \ = by dehnition, Eq. (E136) can be written as 

pi, = (1 - e)^\Ek){Ek\ +[!-(!- e)^] p^. (E137) 

with Ptp orthogonal to \Ek). From Eq. (E134) it follows, 

Nq 

Weye = iVW'ext = ^ [liH) + 0 (/(<?))] , (E138) 

Ph 

where in the last line we have used Eq. (E109). We now set 

N = N{g) = ^— (E139) 

7(k) 9 

for all p > 0 satisfying the constraint N{g) G N+. For any positive constant > 0, one can always consider the values of 
> g > 0 so that N(g) is large. This constraint imposes g = /3hW/(j(K) N), where N has to be an integer. Therefore, 
g now belongs to a subset of the positive real line, rather than the positive real line itself as previously. However, since g 
monotonically decreases to zero as N increases to inhnity, we can still take the limit p —0+ as before. Thus achieving 


Wcyc = w + e (/(p)) . (E140) 

Since limg_j,Q+ /(p) = 0, we conclude part 2) of Corollary 2. For the entropy of the hnal state of the battery we have 

5(p^) = NS ((1 - e) \Ek){Ek\ + e\E,){Ef) = ^ {1 - e)Hl - e) + sine ^ ^ 

As stated above the efficiency is given by Eq. (E108), and thus we can always choose R G (0,1), and p (recall e —?> O'*" as 
p —O'*') such that 1) in Corollary 2 is satished. Furthermore, recall from the proof of Lemma 12 that 


, e\ns 
lim - = 0, 

3^0+ p 


(E142) 


for all K G (0,1). Thus, from Eq. (E141) we conclude that 3) in Corollary 2. We will now prove part 4) of the Corollary. From 
Eq. (E137) and part 4) of the Corollary, we can identify r = (1 — e)^. We thus study the limit 


lim (1 - e)^ 

g^ 0 +' 


f lim (1 

\c)-> 0 + 



PhWhiK) 


( lim 

S-*' 0 + 



PhW/^(K.) 


= 1 , 


(E143) 


where going to the last line, we have used that fact that Eq. (E142) implies that e/p—>^0asp—?'0+. We thus conclude part 4) 
of the corollary. □ 

Thus by choosing 6 > 0 sufficiently small in Corollary 2, we can extract any hnite amount of work with an arbitrarily small 
entropy contribution with an efficiency arbitrarily close to the Carnot efficiency as long as H < 1. 
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F. Extensions to the setup 

Arguably, one may think that the inability to always achieve the Carnot efficiency in the nano regime is due to some subtlety 
of our setup (even though we have shown that according to the standard free energy one can always achieve the Carnot efficiency 
with our setup). For such reasons, in the next few sections we show that even under more general conditions than those laid out 
in Section A, one still cannot achieve the Carnot efficiency when fl > 1. 

In Section F 1, we show that allowing for correlations between the final state of the battery and cold bath (and/or the finite 
dimensional machine) does not allow us to achieve the Carnot efficiency. The main result is Theorem 3. 

In Section F2, we show that allowing for the battery to be any state with trace distance e from \Ek){Ek\y^ cannot allow us to 
achieve the Carnot efficiency when 17 > 1. This shows that whenever we are unable to achieve the Carnot efficiency, it is not a 
artificial defect from an overly specified battery model. The main result is Theorem 4. 


1. Final correlations between battery, cold bath, and machine 

In Section E 2 a, we stated that the final state of the heat engine after tracing out the hot bath was of the form 

trHot(PcoldHotMw) = Pcold C) Pm C) Pw (FI) 

where p^ = e\Ej){Ej |w + (1 — £)\Ek){Ek\w, i-e. the final state of the charged battery was a tensor product with the cold bath. 
We also demanded that the heat engine is cyclic i.e. that pj^ = p^. In this section, we show that if one allows for the final state 
of the battery, cold bath and machine to become correlated^, one still cannot achieve the Carnot efficiency when 17 > 1. That is 
to say, in this section we allow the final state to be 

tl'Hot(PcoldHotMw) = PColdMW (F2) 

with only two natural constrains, namely that our heat engine actually extracts work, i.e. that 

p^ = e\Ej){Ej\\f^ + (1 ~ £)\Ek){Ek\w, (F3) 


as before, and also that the heat engine is still cyclic, i.e. 


Pm = Pm- (F4) 

Throughout this section, (unless stated otherwise) we will write PcoidMW refer to any generic tripartite quantum state on the 
cold bath, machine and battery satisfying Eqs. (E3) and (E4). 

• In Section E 1 a, we first define the generalized efficiency where one is allowed to consider correlated final states. We 
see that although this may potentially affect the amount of extractable work Wext, the amount of heat change in the bath 
remains the same, by making use of energy conservation and the fact that the global Hamiltonian iTcoidHotMW does not 
contain interaction terms between subsystems. 

• In Section E 1 b, we make use of the generalized second law when a = 1 (which is also the macroscopic second law), in 
order to show that final correlations still do not allow the surpassing of Carnot efficiency. This can be proven by noting 
that the von Neumann entropy is subadditive, and the result is summarized in Lemma 19. A proof sketch can be found in 
the beginning of Section E 1 b. 

• In Section E 1 c, we turn to the case where f7 > 1, where without final correlations it is shown in Theorem 2 that Carnot 
efficiency cannot be achieved. 


a. Defining the generalized efficiency 

Recall that before (see Section C2), we have shown in Eq. (C17) that if the following assumptions hold; 
(i) the final reduced state of the battery p^ is fixed by Eq. (A6), 


^ Recall that the final state of the cold bath, machine and battery are already allowed to become correlated with the hot bath 
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(ii) the state of the machine is preserved, 

(iii) the final state is of tensor product form, i.e. PcoidMW = /°Coid ® Pm ® Pw’ 

then the efficiency for a particular transformation PcoidHotMW PcoidHotMW simplifies to being only an explicit function of Pcoid 
instead of the global final state. This simplified expression of the efficiency in Eq. (C17) is then used to evaluate, for example, 
P™^'^(Pcoid) ill (C5). Since we now drop Assumption (iii) for the final state being uncorrelated, the efficiency and the work 
extracted Wext will now depend on the tripartite final state p^oidMW instead. 

Therefore, let us first write a generalized expression for the maximum efficiency corresponding to a transition PcoWHotMW 
PcoidHotMW ''i^ unitary operator U (t) in this generalised setting: 

P''™(PColdMw) •= ®np ?7(pcoi(j, IText) S.t. trHot[t^(f)PColdHotMwC^(i)^] = PColdMWt 
Wex, 

[Uit),H]=0, (F6) 

= e\Ej){Ej\^N + (1 — e)\Ek){Ek\w, (F7) 

Pii = pSi- (F8) 


See Fig (1) in main text for a definition of the other quantities appearing in Eq. (F5). Recall that the definition of rj is given by 
ry = /S.EI as in Eq. (Cl). In Section C 2 we showed that this can be simplified to 

77 = (l-e + AC/IEext)-\ (F9) 

where AC = AC(pcqJ(]). This equation holds under Assumption (i) and (ii), together with the fact that the global Hamiltonian 
does not contain interaction terms between both baths, battery, and machine. Since the derivation of Eq. (F9) does not require 
Assumption (iii), it still holds for a general tripartite final state PcoidMW- However, dropping Assumption (iii) may potentially 
allow for larger values of Wext, and therefore subsequently might affect For this reason we write p''™ = p‘’™(PcoidMw) 
remind ourselves that it is a function of the entire final state PcoidMW- 

We have written p = p(Pcoid> H4xt) to explicitly show the Wext dependency of rj. Although not written explicitly in Eq. (F5), 
we should remember that U (f), iTnot and 77 m are arbitrary, other than satisfying condition (A.4) in Section A. As such, by 
maximizing ry over Wext, these quantities will accommodate their optimal values to maximize p‘'™(PcoidMw)^' Throughout this 
section, we analyze Eq. (F5) only in the case of near perfect work (Def. (2)) since the proof that perfect work is not possible 
(see Lemma 8 ) also applies to Eq. (F5)^. 

For the purpose of our proofs, we need to define a new family of intermediate efficiencies. They provide the maximum 
possible efficiency, when considering only a particular instance a > 0 of the generalized second laws. For any a G [0, oo), let 
us denote 

7a™(PColdMw) = sup p(Pcoldi H4xt) S.t. Ea(T^old ® Pm ® PW) "^ColdMw) ^ -^^(pcoldMWJ "^ColdMw): (FIO) 

W„, 

tr(-HtPcoidHotMw) = tr(77t PcoidHotMW) I (Fll) 

p^ = £|7;,)(7 ;,|w + (1 - e)|£;fe)(7;fc|w, (F12) 

Pm = Pm- (FI 3) 

See Eq. (B 8 ) for definition of Fq,. We denote ry^ = lim ry)’™. The condition Eq. (Fll), is always satisfied when all the second 

O'—>-oo 

laws are satisfied. We add the condition as a constraint here, since we will need it in order to write the efficiency rj in the form 
ofEq. (F9). 


b. Final correlations do not allow the surpassing of Carnot ejficiency 

In this section, we first show that Carnot efficiency cannot be surpassed even when we allow arbitrary final correlations in the 
final state PcoidMW- This can be done in the following steps: 


This is an advantage, since it rules out cases such as when the Hamiltonian does not support a thermal state (e.g. when the corresponding thermal state’s 
partition function diverges). In this section we consider any cold bath Hamiltonian i?cold that satisfies (A. 6 ) in Section A (i.e. finite dimensional). As such it 5 


will always h 
For the sake ( 
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1. Using the definitions of generalized efficiency (allowing correlations) in Eq. (F5) and generalized intermediate efficiencies 
in Eq. (ElO), we prove an inequality between fi‘’™(PcoidMw) ^“"(PcoidMw)’ a > 0. This is done in Lemma 17. 
From this, we also conclude that ?7‘>”(PcoidMw) < vTIpcoUmw)- 

2. On the other hand, we show that for any final state of the cold bath, machine and battery PcoidMW’ generalized inter¬ 

mediate efficiency for a = 1 only increases, if we consider the tensor product of the marginals PcoidMW- other words, 
7 i™(PcoidMw) — 7 i™(Pcoid ® Pw ® Pm)' ‘'^*1 intuitively see why this is true: it comes from the fact that the von 
Neumann entropy is subadditive, therefore the final state Pc^i^j 0 Pw ® Pm contains more entropy than Therefore 

according to the a = 1 second law, one can potentially draw more work by going to the state Pco^ ® Pw ® Pm instead of 
a con-elated state PcoidMW- 

3. Since the argument for Pi™(Pcoi(i <E>p^<S> Pm) is of tensor product form. Assumption (iii) holds as before, and therefore the 
efficiency only depends on the final state of the cold bath Pc^y. This means that Eq. (FIO) for a = 1 reduces to Eq. (C5). 
Lastly, by using Lemma 22, this allows us to further show in Lemma 19 that even by allowing conelations in Pcoi^mw’ 
efficiency p‘'™(PcoidMw) never surpass the Carnot value. 

Firstly, let us fix the following notation: for an i?-partite state PAiA 2 ...Aij> define the uncorrelated counterpart 

R 

PAiA2...Afl := (R)Pa.' (F14) 

i=l 

Comparing paiA 2 ...Ah and paiA 2 ...Ah, one will see that each subsystem has the same reduced state, but the global state is 
different. Another useful thing is to note that if one is given a Hamiltonian which does not contain any interaction terms between 
each subsystem, i.e. 


R 

HaiA2...Ah — 'y ) 1 ai < 8 ) ■ ■ • Hai • ■ ■ Iah) 

i=l 


(F15) 


then we may conclude that 


R R 

tr(iTAiA2...A„PAiA2...A„) = ^tr( 7 TA.PAj = ^tr(iTA,P^) = tr(iTAiA2...AHPAiA2^)' (F 16 ) 

z—1 i—1 

Lemma 17. For all a> 0 and all states PcoidHotMW’ 

P''™(PColdMw) ^ 7a™(PColdMw)! (F17) 

where 77 “'™ and pS™ are defined in Eqs. (F5) and (FIO) respectively. 

Proof. For every a > 0, Eq. ® ph® Pw^^-coidMw) > Fc{Pco\dM'N^ in Eq. (FIO) is a necessary condition 

for the transformation PcoidMW PcoidMW *^0 occur under an energy preserving unitary with the aid of a catalyst [5]. Energy 
preserving unitaries also preserve the average energy and thus the Eq. in Eq. (FIO) is also 

a necessary condition. If a unitary U{t) satisfies the conditions in Eq. (F5), then by the second laws it satisfies Eq. (FIO) for any 
particular a > 0. As a consequence of these observations, the set of allowed unitaries U(t) in Eq. (F5) is a subset of allowed 
unitaries facilitating the catalytic thermal operation which transforms PcqIjjmw to Pc^i^mw Eq. (FIO). □ 


Lemma 18. For any final state PcoidMW> consider the quantity Pi™(PcoidMw) defined in Eq. (FIO). Consider the optimization 
problem 


®(PColdMw) ■— SUpp(pcold! W4xt) S.t. (tcoU 0 ® Pw: 'fiC^aldMw) — 'Fl(pcoldMWi '^ColdMw)) 




tr (Ht PcoldHotMW ) PcoldHotMW) ) 

Pw = £|-E'j)(£’ilw + (1 — £)\Ek) {Ek\sN 1 
Pm = Pm- 


(E18) 

(E19) 

(F20) 

(E21) 


Then, t]\ (PcoidMw) ~ “(PcoidMw)- 
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Proof. We begin by noting that the free energy Fi can be written as 

F,{p,T'^) = tx{Hp)-^^S(,p), (F22) 

where {H)p := tr[fFp], and S{p) = —tr(plnp) is the von Neumann entropy, while is the thermal state at inverse temperature 
Ph for the Hamiltonian H. Also, let us recall that Wext = — Ej^ > 0 where Ej^ is a constant. 

Next, we consider the free energies Fi ® Pm ® '^CoWmw) ^"d Fi (PcoidMW) '^CoWmw) respectively, and how they relate 
to IFext- First of all, note that the quantity Fi ® Pm ® Pwi rcoidMw) simply a constant that does not depend on Wext- This 
is because 


El ('Tcoid ® Pm ® Pw) '^CoidMw) — El ( 


= Ei{ 
= Ei{ 


T-O 

"Coldt 

T-0 

"Coldt 

T-0 

"Coldt 


'^Cold) + T’l('rMi 

'^Coid) + T’i('rMi 
'^Coid) + T’i('rMi 


Tm) + Ei{t^,t^) 

4)+tr(HwP^)-/3,-'5(p^) 

r^) + Ef, 


(F23) 

(F24) 

(F25) 


where the first two terms do not depend on the battery Hamiltonian at all, while in the last equality we have made use of the fact 
that p^ = \Ej){Ej On the other hand. 


T"! ( PcoldM W t r'coldMW ) — ^r (-ffcoldMW PColdMW ) i^/t ^ (PColdM W ) 


= tr(iTcoidpLd) + tr(iFMPM) + tr(fFwPw) “ Ph ^'S'(PcoidMw) 
= tr(iTcoldPcold) + tr(iFMPM) ~ Ph ^‘^'(PcoldMw) + ^EJ + (1 


e)E'^- 


(F26) 

(F27) 

(F28) 


Note that again, tr(fFcoidPcoid) ^^{EmPm) do riot depend on the battery Hamiltonian and therefore do not depend on E'^. 
Similarly, *S'(PcoidMw) depends only on the eigenvalues of the state, and is independent of E'^. Since e S [0,1), we may 
conclude the following: -^(PcoidMW;’^CoidMw) continuous function that strictly increases w.r.t. E'^, and therefore it also 
strictly increases w. r. t. Wext- 

To prove this lemma, it suffices to show that the supremum over Wext in Eq. (FIO) for a — 1 has to be achieved when 
T’i('rcoid ® Pm ® Pw^r-coWMw) = T’i(PcoidMW> rcoWMw)- We prove this by contradiction. Suppose that Wext achieves the 
supremum for p)*”, and for this value of Wext, -Fi('rcoid ® Pm ® Pw>rcoidMw) > T'i(pcoidMW: r-coidMw)- Since we know that 
E'(pLidMW’'^CpidMw) strictly increases w.r.t. Wext, there must exist an W^j > Wext such that Fi(r^eid 0 Pm ® PwttcoidMw) ^ 
T"i(PcoidMW: rcoidMw)' Furthermore, since by Eq. (F9) we know that the efficiency is monotonically increasing w.r.t. Wxt as 
well, it follows that W^t achieves a higher value of efficiency compared to Wext while satisfying the required constraints at the 
same time. This is a contradiction, and therefore we conclude that the optimization for p)'™ can be simplified to a(PcoidMw)’ 
where the constraint on Fi holds with equality. □ 


Lemma 19. For any final state PcoidHotMW> Hamiltonian of the form in Eq. (Al), then for perfect or near perfect work 

extraction (see Defs. 1 and 2 ), we have 

rf^{PcMMw) < Pr(PMdMw) < pr ( p^idMw ) = P™" (p^id) < 
with equality in (2) (j^PcoidMW = PcoidMW- TPe quantities p)*™ and are defined in Eq. (FIO) and Eq. (C5) respectively. 


Proof. Note that inequality (1) is a direct consequence of Lemma 17, while inequality (4) holds because of Lemma 6. It remains 
to prove inequalities (2) and (3). 

Proof of inequality (2): Using the definition in Eq. (LIO) together with Lemma 18, let us compare the quantities 

Pi (PColdMw) = SUpp(pcold) ^ext) S.t. -Fl (tcoW ® Pm ® PW) "^ColdMw) = T’l (PColdMW)'^ColdMw)) (L30) 

Wex, 

tr (Ht PcoidHotMW) = tr (PcoWHoimw ), (F31) 

p^ = e|F,)(F,|w + (1 - e)|Ffc)(Ffc|w, (L32) 

Pm = Pm, (F33) 
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and 


^1 (PColdMw) — ^(Pcold) ^ext) S.t. Fi 0 Pm ® PW) "^ColdMw) — ^1 (Pc 


Wex 


ColdMW; 'ColdMW 


), 


tr(-^iPColdHotMw) — tr(-ffiPcoldHotMw)) 

Pw = £|^i)(-E’ilw + (1 ~ s)\Ek){Ek\yN^ 


Pm — Pm- 


We first make the following observations: 


(F34) 

(F35) 

(F36) 

(F37) 


• By our definition of PcoidMW’ *^Fat pJoid = Pcoid- Therefore, the term AC in Eq. (F9) which is only a function 

of the reduced state on the cold bath is the same for both efficiencies in Eq. (E30) and Eq. (E34). Therefore, to compare 
the efficiencies, we need only to compare the value of W^xi that satisfies the free energy constraint in both optimization 
problems. 


• In [30] (pg. 395) it has been proven that the von Neumann entropy is subadditive 

S{pab) < S (p^) , (E38) 

with equality iff = Pab- Eurthermore, since TT^coWmw does not contain interaction terms, as we have demonstrated 
earlier in Eq. (E16), 

tr(f7coldMWPColdMw) = tr(-ffcoidMw PcoidMw ) ■ (E39) 

Thus, by Eq. (E22) we conclude that 

F’l( PcoldMw ) — -Fl(pcoldMw)) (E40) 

with equality iff p^^^mw = PcoidMW- 


• Eor any final state PcoidMW where p^ = e|£’^)(Cj|w + (1 — e)|Cfe)(i?fc|w, we have seen in the proof of Lemma 18 that 
-Fi(PcoidMWJ '^ColdMW) ^ continuous function that strictly increases with Wext- 


With these three observations we can now prove inequality (2). Note that when PcoidMW = PcoWmw’ equality holds trivially. 
Therefore, let us consider the case where Pcommw PcoidMW- Suppose lEext achieves the supremum in Pi™(PcoidMw)’ ^^'d for 
such a value of VFext, 


-Fl('rcold Pm Pwi 7'ColdMw) — -Fl(pcoldMWJ '^ColdMw) > F'i(pcoidMW> "^ColdMw)- 


(F41) 


We note also that since F’i(PcoidMWi 'h^oidMw) strictly increases with Wext, and therefore there exists some Wg'^t > l^ext such that 
El (Tcoid ®Pu® Pw. '^CoidMw) = El ( PcoidMW : '^CoidMw)' Therefore, is a feasible solution for Eq. (E34), i.e. it satisfies the 
constraints in the optimization problem. In conclusion, we have 


Pi (PColdMw) — 


1 — e - 


AC 


n -1 




< 


1 - e + 


AC 


n -1 


^ Pi (PColdMw)- 


(F42) 


Proof of equality (3): Consider the quantity Pi™(PcoidMw)- Since the state PcoidMW takes on a product structure form between 
all the subsystems now. Assumption (iii) in the beginning of Section E 1 a holds again. By Eqns. (E32) and (E33), we know that 
Assumptions (i) and (ii) also hold. Therefore, we know that under these assumptions the efficiency does not depend anymore on 
the global state PcoidMW’ Fut only Pcow Again comparing the conditions of p™^‘^(pcoid) ^ttd Pi™(PcoidMw)’ we see that they are 
exactly the same quantity. □ 


Therefore, Lemma 19 tells us that correlations between the hnal states of the cold bath, machine and battery cannot allow you 
to achieve an efficiency greater than the Carnot efficiency. 
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c. Achievability of Carnot efficiency still depends on more than temperature 


Earlier in Section F 1 b, we proved in Lemma 19 that Carnot efficiency gives an upper bound to the efficiency of any arbitrary 
final state PcoWMW- '^^is section, we want to prove that when > 1 holds, one cannot achieve the Carnot efficiency even when 
allowing correlations between the final states of the battery and the cold bath. This can be done in the following steps: 

• According to Lemma 19, Carnot efficiency can be attained only when all the inequalities in Eq. (F29) are satisfied with 
equalities. We use this to prove in Lemma 20 that in order to achieve the Carnot efficiency, we may only consider the limit 
where correlations in the final state vanish. Not only so, the magnitude of these correlations also have to vanish quickly 
enough in order for Carnot efficiency to be achieved. In particular, we define a parameter k which quantifies the amount of 
correlations, and show that k has to vanish faster than the quasi-static parameter g, in order to achieve the Carnot efficiency 
VC- 


• Next, in Lemma 21, we show that if the parameter k vanishes faster than the quasi-static parameter g, then whenever 
n > 1, one can derive an upper bound for the intermediate efficiency Poo’(PcoidMw) which considers the amount of work 
extractable by invoking only the generalized second law of a —oo. Combining Lemma 20 and Lemma 21, we conclude 
in Corollary 3 that when > 1, p'’™ < < rjc is strictly upper bounded away from the Carnot efficiency. 

Before we begin, let us note that by definition, the initial state Pcoidw i^ diagonal in its energy eigenbasis. Furthermore, the 
state PcoidMW i® of form Pcoid ® Pm ® Pw- Since w.l.o.g. we can assume that Hm is proportional to the identity (or called the 
trivial Hamiltonian, see [5]), can always be written as a diagonal state in an energy eigenbasis of its Hamiltonian. Therefore 
the state PcoidMW i^ always diagonal in the energy eigenbasis of the Hamiltonian ITcoidMW := f^Coid + Hm -b Since catalytic 
thermal operations cannot create coherences [5], PcoidMW f’o also diagonal in the energy eigenbasis of ITcoidMW- 

We observe that any PcoidMW ^an always be written as 


PColdMW — (1 ^ )PColdMW + ^ PColdMWi 


(F43) 


where k* = min{/c S [0, l]|PcoidMW = (1 “ ^)PcoidMW + kQ,Q > 0}. This means that PcoidMW can be written as a convex 
combination of two states: one being PcoidMW’ other p^oUMW containing all other correlations. Note that such a k* always 

exists, in particular, fc = 1 is always a feasible solution. 

We now dehne a particular parametrization of the final states, 

PColdMw(^: P&ldMWi PoddMw) ■= (1 ~ ^)PCoWMW + ^PoddMWi ^ G [Oj (F44) 


where the following holds: 


/•\ no corr _ 1 

PColdMW — PColdMW^ 


con' / no con 

PColdMW ^ PColdMW^ 

(hi) = + = 


(F45) 

(F46) 

(F47) 


Since in our heat engine, the initial state has no coherences, it suffices to consider PcoidMW which is diagonal in the energy 
eigenbasis. This implies that PcoidMW = PcoidMW diagonal in the energy eigenbasis, and therefore the same holds for 

PoddMW (F44). All correlations between the individual systems of cold bath, machine and battery are contained only 

in PcokiMW- Therefore, PcoidMw(’i ’) parametrizes every possible quantum state on 'HcoidMW which is diagonal in the global 
energy eigenbasis and that returns the machine locally to its initial state after one cycle of the heat engine. In Eq. (F47), Pm is 
the final state of the machine, since the heat engine is cyclic, recall from Section A that we require Pm = Pm- 

Lemma 20. For every family of states PcoidMw(^! PcoidMW! PqJwmw) parametrized by k, (see Eqs. (F44)-(F47)), if the quantum 
efficiency 77 ^™ defined in Eq. (FIO) achieves the Carnot efficiency 

Vi (PcoidMw) = 1 ~ (F48) 


then the following conditions are satisfied: 

1) The state PcoidMW the final state of a quasi-static heat engine (see Def 3) 


PcoidMW =t(p)®Pm(p)®Pw withg^ 0 +. 


(F 49 ) 
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2) The correlations must vanish sufficiently quickly. That is to say, the parameter k in Eq. (F44) vanishes more quickly 
compared to g, i.e. 


k 

lim - = 0. (F50) 

3^0+ g 


Proof. Firstly, suppose that Carnot efficiency is achieved, i.e. ??‘*™(PcoidMw) = 1 ~ Then according to Lemma 19, all 
inequalities in Eq. (F29) should be satished with equality, in particular inequality (4). We have established in Lemma 6 that this 
equality is only achieved in the quasi-static limit, i.e. Pcoid ~ ''Coid(ff) where p —O’*". This implies Condition 1) in the statement 
of the lemma. 

The proof for Condition 2) consists of calculating Wext for a = 1 in Eq. (ElO) to leading order in g and k. This ILext quantity 
can be later used to evaluate ryj™. We will show that we can write the expression for pj™ into two terms: one term describes the 
efficiency when there are no hnal correlations, and the other term is a strictly negative contribution which must vanish in order 
to achieve the Carnot efficiency. This latter constraint will give us Eq. (E50). 

Let us denote as the value of battery energy gap Wext = E'k ~ that solves the equation 


T’l('^Cold C) Pm( 5) ® Pw) tColdMw) — T’l(pcoldMw(^) POtUlMW) PoddMw)) '^ColdMw)^ 


(F51) 


while Wext as the value that solves the case where k = 0, i.e. 

T'l('tCold ® Puid) ® Pwt fColdMw) = T’l(PcoldMWi ’^ColdMw)' (F52) 

Since PcoidMW = PcoidMW = Pcoid® Pm( 5') contains no correlations, Wext was given by Eq. (D51). According to Lemma 18, 
we know that Wext and Wext are the values of Wext which solve sup ??i™(PcoidMW) fFext) and sup ??i™(PcoidMWi fFext) respectively. 

Wext Wex, - 

Solving Eq. (E52) for Wext with the aid of Eq. (F22), we hnd 

FFe'xt = fFext - X, (E53) 


where Wext is the solution to Eq. (F52) when k = 0, given by Eq. (D51), while 

F ■“ ^ ['^'(P&rdMw) ~ ^ (PColdMw(^W&ldMWi PoddMw))] • 

Let us hrst note some properties of %, which we will later use: 

• Since S{-) is subadditive, due to the parametrization of PcoidMw(’t’) i^^ Fq. (F44), we have 


(F54) 


X > 0 


with equality iff PcoidMW ~ PcoidMW i-C- iff ^ = 0. Therefore, we may conclude that > 1. 
• We have that 


if and only if 


d 

dk 


X(fcj P&tldMW) PaJidMw) 


= 0 

k—ko 


PColdMw(^0i PcoldMWt POtkiMw) 


IcoldMw/-^- 


(F55) 


(F56) 


(F57) 


Eqs. (F56) and (E57) are direct consequences of the observations: 

1) Entropy is strictly concave, i.e. S (PcoidMw(^W&idMwWa)'idMw)) strictly concave in A: G [0,1]. Therefore, by 
Eq. (F57) X is strictly convex in fc G [0,1]. When the hrst derivative of the convex function ^ — 0, this must be the 
global minimum [? ]. 

2) However, we know that the entropy is uniquely maximized (and therefore x is uniquely minimized) for the maximally 
mixed state. 


® We denote p^(g) because for different values of g, we are allowed to choose different initial machine states, as long as pjj(3) = p^{g)- 
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Returning to evaluate the efficiency, we may use Eq. (F9) to calculate the inverse efficiency. 


[7r(PColdMw)]-' = l-e+^^^ 

^^ext 

(F58) 

_ 1 Wext 

Wext W'xt 

(F59) 

- 1 AC'(pcoid) 

(F60) 

Wext 


The last term in Eq. (E59) is non-negative because we know the terms AC{pQgi^), Wext and are all non-negative. 
With Condition 1), we now know that 


Wext P ’ 


in the quasi-static limit, and therefore a necessary condition to achieve the Carnot efficiency is that lim„_i.o = 1 also in the 

^ext 

quasi-static limit. Using the relation = VEext + X, we have the requirement that 


lim 

9->0+ 


X(fei PcoidMw(g)> POildMw) 

Wext (PcoldMW (5 ) ) 


= 0 . 


(E62) 


First, let us observe that Wext (PcoidMW( 5 )) = Wext(/3c-5) given by Eq. (D37). The leading order term of Wext( /3c -ff) = ©( 5 ) 
as g —>■ 0+. Therefore, in order to satisfy Eq. (F62), we must firstly have limg_i.o X = 0- From Eqs. (F44), (F54), this implies 
that we need A: —?► 0 for all PcoWMW 

Since the numerator and denominator of Eq. (F62) both go to zero, by L’Hospital rule, to evaluate the limit we need to take 
the derivative of both terms w.r.t. g. Therefore, we expand x to first order in k and g. From Eq. (F54) it follows 


xik, P&ldMW ( 5)1 POikiMW) “^^^’ PColdMW(0); POikiMW) 


k—0 


^X(0i PMdMw(ff)j PoddMw) 


g + o{gk) + oik"^) + oig"^) 


9=0 


— ^ X ( ; PColdMW ( 0 ) ! POikiMW ) 


k + o{gk) + o{k'^) + o(p^). 


fe =0 


(F63) 

(F64) 


The term j|x(0, /0^°idMw(ff), PcoMmw) 


0 since when 


9=0 

(F55) holds, it must be that ^x{k, PcoidMw(O): PoddMw) 


k = 0, X will be constant for all g. Next, we note that since Eqs. 
> 0. Furthermore, from Eq. (F56), we have that 


k=0 


^ 0, (F65) 

k^O 

for all pgJdMW since by definition p^oidMw(0. PcoidMw(0)> PSJdMw) 7^ 1 CoidMw/3V. We can infer that *s not maximally 

mixed from a few observations, for example: this is true because we have required that the reduced state on the battery is not 
maximally mixed since we consider near perfect work extraction. 

Thus, taking into account ~ 9) = ©(5)> Eq. (F62) implies Eq. (F50). □ 


^ X ( ^ I PColdMW ( 0 ) I PoddMW ) 


By now, we have established a constraint on how quickly the correlations have to vanish w.r.t. g, for the possibility of achieving 
Carnot efficiency. In the next Lemma 21, we will show that the constraints given by Eq. (F50) can be used to derive an upper 
bound for 
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Lemma 21. IfEqs. (F49) and (F50) are satisfied, then the quantity can be upper bounded by 

^oo (PColdMW ( ^ J P&WMW (5 ) ) PaJklMW ) ) 


< 


1 


Ph 7(1) 


1 -1 


Pc - Ph 7(00) 


+ Q{f{g)) + e[k/g) + Q{g) + eie) 


(F66) 

(F67) 


with liing^o+ fid) = 0. 


Proof. The main idea of our proof is as follows: we show that if Eqns. (F49) and (F50) hold, then we can upper bound Wgxt 
while considering only the Fao condition. This bound differs from the value given when no correlations are present by only a 
small amount. Substituting this into the expression for rf^, we obtain Eq. (F67). 

Let us begin by analyzing the difference in eigenvalues of the states PcoidMW /^CoidMW Recall that 


PColdMwC^J P&ldMW! PoddMw) 


(1 — fc)PcoldMW + ^POddMW 


(E68) 


where PcoidMWi /^SkiMW diagonal in the energy eigenbasis. Since PcoidMW ^ mixture of two energy-diagonal states, it 

is also diagonal. Let us denote its eigenvalues as [PcoidMw]*- 

As for PcoidMW’ ^9*'- (F49) gives the explicit form of the state. 


PColdMW — Pcold ® Pm ® Pw — '''(p) ® Pm(p) ® Pw- 


(E69) 


Let us denote its eigenvalues as [PcoidMw]*- 

We hrst observe two properties involving trace distance 


(P.i) Consider two states ci, 172 diagonal in the same eigenbasis. Then if p = (1 — k)ai + k(j 2 for some k S [0,1], then one 
can conclude that the distance 


d{p,cri) < k. 


(L70) 


(P.ii) Lor any two states p, a diagonal in the same basis, with eigenvalues pt, qt, if their trace distance 

dip,cr) = ^\\p-(^\\i < £, (F71) 

then this implies that their eigenvalues cannot differ by more than e, i.e. \/i, \pi — qi\ < e. By using this fact, we may hrst 
calculate the trace distance between PcoidMW PcoidMW’ *^hen bound the difference of their eigenvalues. 


We hnd that 


^(PColdMWi PColdMw ) — '^(PColdMWi PColdMw) + ^(PColdMWi PColdMw) 

<A + rf(pcoir,Pcoid)+rf(P' 

< Ak. 


no corr 


,PM)+rf(Pw“",Pw) 


(P72) 

(P73) 

(P74) 


The hrst inequality is a triangle inequality that holds for all states. The second inequality holds because of (P.i), and because trace 
distance is subadditive under tensor product (note that both PcoidMW PcoWmw tensor product states). The third inequality 

holds because we know (i(PcoidMW) PcoidMw) — ^ '^^at trace distance decreases under partial trace. By (P.ii), Eq. (L74) tells 

us that \/i. 


[PColdMw]* — [PColdMw]* + o{k). 


(L75) 
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With Eq.(F75), we may relate the Foo quantities for the states PcoidMW PcoidMW- From Eq. (B 8 ), we have 


Foo (PcoldMw(^i PColdMw(ff)) PoddMw): r'coldMw) 

[PColdMw]* 


= In max 

i 

= In max 


[PColdMwli 


Fo(fc), 

= Foo (Tig) 0 Pm(5) ® Pwj 'rcoidMw) + o{k), 


n 

„0 


where we used Eq. (B 8 ) in the last line. 

The next step is to evaluate the restriction on Wext that satisfies 

Foo('rcoid ® Pm ® Pwj rcoWMw) ^ ^00 (PColdMw(F, PcoldMw(p)) PoddMw)) r'coldMw) 

= F'oo (t(p) 0 Pm(p) ® Pwi rcoidMw) + o{k), 


(E76) 

(E77) 

(E78) 

(E79) 


(E80) 

(E81) 


for Wext up to order o{k). Taking into account the additivity of F^o under tensor product, we can rearrange Eq. (E81) to provide 
an upper bound on Wext, 


Wext < ^ [ 7 ( 00 ) + <d{f{g)) + o{k/g)], 


(E82) 


where lim f(g) = 0, 7 ( 00 ) is given by Eq. (El 13). The bound in (F82) is achievable since the F^o conditions imposed by Eq. 

3^0+ 

(E80) are achievable with equality. 

Lastly, by using the expression for efficiency in Eqs. (F9), and substituting Wext from Eq. (F82) (with equality for the 
maximum possible Wext) followed by AC from Eq. (El 16), we have 


sup p(Pcold> Wext) 
Wm >0 


sup I 1 — £ + 

Wext>0 


AC Y 

Ph 7 ( 1 ) 

(/3c - Ph) 7(oo)_ 


1 -1 


+ 0 (/(p)) + oik/g) + 0 (p) + 0 (£). 


(F83) 

(F84) 


Hence using Eqs. (FIO), (F84) we find Eq. (F67). Note that in Eq (F67) we have an inequality, this is due to the fact that in 
the optimisation problem Eq. (FIO), there is an additional constraint (namely mean energy conservation) which is not taken into 
account in the derivation of Eq. (F84). □ 


Finally, the above lemmas allow us to conclude that allowing further correlations in the final state cannot allow us to achieve 
the Carnot efficiency when H > 1. 

Theorem 3. [Correlations do not improve efficiency] Suppose that H > 1. Parametrizing the final state of the heat engine by 
Eq. (F44)-(F47), the quantum efficiency ffi”' defined in Eq. (F5) is strictly upper bounded by the Carnot efficiency, 

sup p'’™ (PcoldMw(^) P&ldMW: PQddMw)) ^ ~ IT' (F85) 

fee[0.1], pSiTmw Fc 

Proof. From Lemma 17, we have that both < pY' und < rj^ hold. Thus a necessary condition to achieve the Carnot 
efficiency for a particular Pcommw’ that both p)'” and p^ are equal to or greater than the Carnot efficiency. 

Lemma 20 proves that Eqs. (F49) and (F50) are necessary conditions for p)*™ to achieve the Carnot efficiency. However, when 
Eqs. (F49), (F50) are satisfied, then Lemma 21 provides an upper bound on the efficiency p^ in Eq. (F67). 

Now, suppose H > 1. Since it is shown in Eq. (El 14) that 7 ( 1 ) /^{oo) = H, plugging this into the leading term appearing in 
Eq. (F67) 



Ph 7 ( 1 ) ' 
(/3c - Ph) 7 ( 00 ). 


-1 


(F86) 


we have that the quantity p§S* (and therefore also p*!™) is strictly less than the Carnot efficiency 1 — Pn/ Pc- 


□ 













57 


2. A more general final battery state 

For the simplicity of our analysis, we have assumed that the battery is left in the specific final state described in Eq. (A5), i.e. 
an amount of work Wext = Ek — Ej is extracted, except with failure probability e that the battery remains in the initial state 
\Ej){Ej\^. In this section, we show that this is a simplification which can be removed in general, i.e. the final battery state is 
allowed to be any state within the e-ball of |i7fc)(^'felw particular, our result that the Carnot Efficiency cannot be achieved 
when n > 1 still holds. 

In Lemma 22, we show that for any final state of the cold bath Pcoid’ allowing a more general final battery state does not affect 
the amount of work bounded by the E^o condition. We then use this to prove in Theorem 4 that when 12 > 1, Carnot cannot be 
achieved even if we allow a more general battery final state. 

Lemma 22. For any given PcoidJ Pcoid’ Pw — \Ej){Ej\y^, consider the maximum := Ek^ — Ej such that Pcoid ^ Pw 
Pcoid ® Pw allowed by the non-increasing Foo condition (Eq. (B7)j i.e. 

^oo(Pcoldlkcold) + > £’oo(Pcoldlkcold) + ^oo (Pw 11) > (F87) 


with 


Pw — (1 l^fciX^feilw E e \Ej){Ej\.^ . 


(E88) 


On the other hand, consider any battery final state 

pl = {l- e) \Ek,){Ek ,Iw + (F89) 

where is an energy-diagonal state orthogonal to \E}^f){Ekf\^ which may depend on e, i.e. = '^^Pi |i?i)(i?i|w 

Pk^ — 0 and J2iPi — 1- Define W'^ := Ek 2 — Ej such that p^oid Pw ^ Pcoid ® Pw allowed by the non-increasing F^o 
condition, i.e. 


77oo(Pco1(i||TqJ‘[jj) -f i2oo(Pwli'^W^) — F7oo(Pcoldll''’Cold) + F)oo (Pw I i'’w*') ' (E90) 

Then for allO < e < i = [l -f we have = Wf^. 

Proof. Firstly, note that any energy-diagonal state with trace distance d{p'^, \Ek 2 ){Ek 2 Iw) = ^ written in the form of 

Eq. (F89). Rearranging the terms in Eq. (E87), 

£'oo(Pwlkw'') < F7oo(Pwlkw'') + F7oo(Pcoldlkcold) - F?oo(Pcoldlkcold) =: (^^l) 

One can use the definition of Dao in Eq. (B12) to expand the L.H.S. of Eq. (E91), obtaining 

logmax{(l — < A — logZ^''. (E92) 

We know that since near perfect work is extracted, e is arbitrarily small. This implies that for e small enough, max{(l — 

£)gP)i-Efci ^ I _ . 

Similarly, one can evaluate Eq. (E87) to obtain 

logmax{(l - , {epie^^^'}i^k 2 } < A- \ogZ^^. (E93) 

Note that the maximization in Eq. (F93) only picks out the maximum value. In particular, denoting Smax to be the largest energy 
eigenvalue of the battery, then whenever 


(1 — 


or equivalently 


e < 


^ _j_ gP/i(.t?max 



(E94) 


(F95) 


then max{(l — e)e^'*^'‘ 2 , {epie^^^'}i^k 2 } = (1 ~ . In other words, as long as e is upper bounded by Eq. (E95), we 

know which terms attains the maximization in Eq. (E92). However, we also want an upper bound that is independent of any 
limit involving the final state PcoidMW "'o wish to take, or any amount of work extracted (and therefore, we want the bound to be 
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independent of E^^). As such, let us constract the following upper bound e < i where. 


e := inf 

-Efc 2 

W ^>0 


1 + e 


,/3)i(-Ema,x —-Bfc,) 


n -1 


o/3h(-Emax —Sj) 


n -1 


(F96) 


Now, we see that Ek^ and Ek 2 correspond to the solutions for Eq. (F92) and Eq. (E93), which for e < e reduce to exactly the 


same equation. Therefore, = E^^ and hence = W^. □ 

We will use Lemma 22 to prove Theorem 4. But before we proceed, let us fix some notation; we define the efficiency as a 
function of a > 0 : 

vi{ploid) = sup ?7(Pcoid) subject to E„0 (Pw ® PcoXi^ '^Coidw)> (F97) 

Ekj —Ej >0 

and tr(iTtp[i„yHotMw) = tr(iTtPco‘idHotMw)- (F98) 

with J = 1, 2 denoting the final battery state We also define an a independent efficiency; 

P'^(Pcoid) = sup r7(pcoid) subject to (F99) 

Ekj — Ej >0 

-P'a(Pw®'’’Coldi '^Coldw) — Fa{p^ 0 PColdJ ''’Coldw) — 9- (FlOO) 

For any a > 0, and any state > P“'(Pcoid) holds. 


We already know that when fl > 1, for any final cold bath state p^oid, the efficiency p^(pcoid) ts strictly less than the Carnot 
efficiency. Theorem 4 shows that this is also true for p^(pcoid)’ he. when allowing a more general battery final state. 

Theorem 4. [General battery states do not improve efficiency] Consider a heat engine with a cold bath consisting of n qubits, 
and consider the case where > 1 (recall the definition of il in Eq. (E106)j. Then for any final cold bath state Pcoid> t/re 
efficiency P^(Pcoid) strictly less than the Carnot efficiency. 

Proof. Firstly, suppose that fl > 1. By Lemma 15 we know that the infimum is obtained at a = c», and by Lemma 16 we know 
that the efficiency for quasi-static heat engine is strictly less than the Carnot value; 

liinpijTp ) = liinp^(Tp ) < PC- (FlOl) 

g-fO ■' g->0 ■' 

On the other hand, we also know from Lemma 6 that p^(pcoid) oan only possibly achieve Carnot efficiency in the quasi-static 
limit. In other words, for all other final states Pcoid know that Carnot efficiency cannot be achieved. Therefore, it suffices to 
see that in the quasi-static limit, 

lim < lim = lim = lim < pc- (F 102 ) 


The second equality is obtained by noting that for any state p^ojd (and therefore for rp^ ); 

1 . AC is the same for both expressions of efficiency p^(Pcoid) P^CPcoid)- 

2. By Lemma 22, for all 0 < £ < [l -f ITi(p^„,d) = IF^(Pcoid)- 


Hence, from Items 1 and 2, one concludes that p^(pcoid) = P^(Pcoid)- ^Fe third equality in Eq. (F102) comes directly from 
Eq. (FlOl). □ 






